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{3JQ' Abstract 

We present an introduction to the canonical quantization of gravity performed in loop 
Cn , quantum gravity, based on lectures held at the 3rd quantum geometry and quantum 

gravity school in Zakopane in 2011. A special feature of this introduction is the inclusion 
of new proposals for coupling matter to gravity that can be used to deparametrize the 
r^ ' theory, thus making its dynamics more tractable. The classical and quantum aspects of 

^S| . these new proposals are explained alongside the standard quantization of vacuum general 

P^ ' relativity in loop quantum gravity. 

o 

1 Introduction 

Loop quantum gravity aims to formulate a quantum theory for Einstein's classical theory 
^^ \ of general relativity. As in the case of quantum mechanics one takes a given classical 

Jh ' theory, in our case general relativity, as the starting point and then tries to develop the 

- - -' corresponding quantum theory. However, when we aim to quantize general relativity, we 

face two aspects that are particular to general relativity. At first, in contrast to the gauge 
theories that have been quantized in the context of the standard model of particle physics, 
in general relativity the metric itself becomes a dynamical object. For this reason we are 
no longer in the situation that we can consider matter put on a fixed background geome- 
try and then describe how matter evolves and interacts with each other on a given fixed 
spacetime. In general relativity matter and geometry interact, and this highly non-linear 
interaction is described by Einstein's equations. Secondly, in general relativity, due to 
the requirement of general covariance, spatial and temporal coordinates do not posses 
any physical significance. This is also refiected in the fact that general relativity can be 
understood as a gauge theory with spacetime diffeomorphisms as the gauge group. Thus, 
observables in general relativity, as these are gauge invariant quantities, are usually more 
difficult to construct than in other gauge theories due to the complexity of the gauge 
group of spacetime diffeomorphisms. In the context of loop quantum gravity (LQG) these 
two points are taking as a guiding principle for formulating the quantum theory. The first 



point mentioned above explains why a non-perturbative and thus background indepen- 
dent quantization is chosen. A more detailed discussion about this motivation can be also 
found in the article of Ashtekar in these proceedings. Although techniques from ordinary 
lattice gauge theory are adopted for LQG, the representation we end up with is rather 
different from the Fock representation used in ordinary quantum field theories and again 
the role of diffeomorphisms enters crucially here. This representation allows to formulate 
a quantum analogue of Einstein's equation, also called Wheeler-deWitt equations, and we 
will discuss those quantum Einstein equations and their derivation in this review. The 
construction of observables and therefore the description of the physically relevant sector 
of loop quantum gravity can be performed in different ways. One possibility is to extract 
the gauge invariant sector of the theory at classical level and then quantize. Another is to 
follow Dirac's idea and derive the gauge invariant part of the theory at the quantum level. 
Also, we can combine those ideas and only partially reduce the system at the classical 
level. In these lectures we will briefly discuss how one proceeds in those cases and what 
are the properties of the resulting quantum theories. One particular aspecty of our intro- 
duction is the inclusion of new proposals for coupling matter to gravity that can be used 
to deparametrize the theory, thus making its dynamics more tractable. We discuss the 
classical and quantum aspects of these new proposals alongside the standard quantization 
of vacuum general relativity in loop quantyum gravity. We hope that this parallel treat- 
ment makes them more acessible and allows a clear comparison to the standard approach. 

Of course we will not be able to present a complete introduction to loop quantum gravity 
but the idea of the article is to provide a rather brief introduction to the techniques used 
in LQG. For instance we will only focus on the canonical part of the theory in this article. 
A path integral approach to LQG in the context of spin foams will be explained in an 
article by Rovelli in these proceedings. We will only briefly discuss the connection this has 
to the canonical approach reviewed here in section|4]on quantum dynamics. Introductions 
to LQG that contain a much more detailed description can be found in Rovelli's [1] and 
Thiemann's book |.2i . A book on LQG that is particularly addressed to undergraduate 
students has been published by Gambini and PuUin [3 . Furthermore, Bojowald's book [4] 
presents an introduction to LQG with a focus on applications of the theory. Other already 
existing lecture notes on loop quantum gravity can for instance be found in [3 |6l [7l |8l |9] . 

We have structured this article into four main parts. The first part in section [5] in- 
troduces the classical framework, that is needed for general relativity in oder to take it 
as a starting point for loop quantum gravity. In section [3] we discuss the quantization of 
general relativity and explain how the kinematics of loop quantum gravity can be formu- 
lated. Afterwards we discuss the quantum dynamics in section |4] and finally summarize 
and conclude in section [5l 

2 Classical theory 

This section deals with the classical setup for loop quantum gravity. We will start with a 
derivation of the Ashtekar variables for general relativity, that are the elementary phase 
space variables used in the classical theory in section 12.11 Afterwards in section 12.21 we 
present a discussion on the dynamics of general relativity with a foucs on those properties 
of the classical dynamics that are relevant for the quantization later on. Finally, we 
present the classical formulation of two models that have been recently introduced for the 
dynamics of loop quantum gravity in section 12.2.11 and 12.2.21 

2.1 Canonical connection variables 

In these lectures we describe a canonical quantization of gravity. Therefore we need 
to start with a canonical description of the classical theory, that is, a description in 



terms of a phase space, canonical variables, and a Hamiltonian. The canonical variables 
are just coordinates in the phase space, so strictly speaking they are irrelevant in the 
classical description. Their change is just a canonical transformation, and hence does 
not change the physics. But in our context, we have to pay close attention to the choice 
of canonical variables, as it affects the quantum description. Canonical transformations 
rarely correspond to unitary maps in the quantum theory: 



Canonical formulation 1 > Quantum theory 1 



canon, trafo 



(1) 



Canonical formulation 2 > Quantum theory 2 



So the canonical variables have to be chosen with care. How to do this? While any choice 
will ultimately have to be justified by the results, one can be guided by other factors. In 
our case, this is the following fact [TUIITTIIT^ : 

There is a formulation of general relativity in which the phase space is is 
precisely that of SU(2) Yang-Mills theory. 

Before we describe this formulation, let us consider the canonical formulation of Maxwell 
theory as a warm-up. The action is 

S[A] = -\ f F,,F^-^d% (2) 

* Jm 

with A the 4-potential and F its curvature. Now we chose an equal-time surface Sj relative 
to some inertial observer time t. The canonical coordinates are the 4-potential A^ and and 
the conjugate momenta P^^ = SS/SAi^^. We find that the spatial momenta are given by 
the electrical field, P" — — _E°. Moreover, P^ — which implies that Aq is non-dynamical 
in our description. We can then rewrite the action as 

S= f dt f d^x - E^Aa -\{E^ + B^) + AoV ■ E. (3) 

From this form, one can conclude the following: 

1. Aa,E^ are the canonical coordinates. 

2. (E^ + B'^)/2 is the Hamiltonian of the system. 

3. Aq is a Lagrange multiplier. It enforces the constraint V -E — 0, which is just Gauss' 
law in the absence of charges. 

Now we will go through the exact same steps for general relativity. To be able to do 
this, one needs to start from a formulation in terms of a connection and a (co-) frame. 
What is a framed It is simply a basis e/ = 6^9^ of the tangent space at each point of 
space-time. Equivalently, it can be viewed as a point-dependent map R — )■ TpM. If there 
is an internal metric, that is, a metric ri on M"*, the frame defines a space-time metric g, 
by declaring the components of the frame to be orthogonal: 

g'"- := e^e:;ry". (4) 

Vice versa, a space-time metric g defines an orthonormal frame, but only up to S0(3,l) 
rotations in the internal space R''. A co- frame is the same as a frame, but with respect to 
the co-tangent bundle, e^ = ei^dx'^. A frame defines a co-frame, and vice versa, by sim- 
ply pulling indices with the corresponding internal and external metric, or, equivalently, 
inverting the matrix of components. 

With that said, we can write the action of general relativity in terms of a co-frame e^ 
and an so(3,l) connection w: 

S[uj,e]^ Sp[uj,e] + S}i[uj,e], (5) 



with [13] 



up J 



KL 



(6) 
(7) 



Here k — IGttG is the couphng constant of gravity, and /3 is an additional parameter, 
the Barbero-Immirzi parameter to be discussed below. The Palatini action Sp was long 
known as a viable action for general relativity. The Hoist action S'h[w, e] does not change 
that: The equations 6S/Suj are equivalent to Z3^"'e = [13], which can be solved for 
id = uj(e). Re- inserting w(e) into the action gives 



S'H[w(e),e] = 0, Sp[oj{e),e] = S'EHL9(e)], 



(8) 



where S'eh is the Einstein-Hilbert action. Thus the above action principle leads to the 
equations of motion of general relativity, irrespective of the value of /3l^ The values 
/3 = ±i are, however, special in the sense that the resulting canonical formulation has 
special properties [lOj . as we will see shortly. 

To go over to the Hamiltonian formulation, we first have to chose a time function 
t on the space-time manifold, whose level surfaces St give a foliation of the space time 
into spatial slices. Additionally we pick a time co-vector field t" with f^dat = 1, and 
decompose it into tangential and normal components with respect to Sj, 



r = Nn" 



N° 



(9) 



where n" is the unit normal, and the shift vector N'' is tangential. Next, we partially fix 
the gauge freedom in the co- frame by linking it with the normal one- form «„, 



e° = "a- 



(10) 



This gauge is called time gauge. Since n^ is time-like, only S0(3) remains as gauge group. 
The covariant fields can now be decomposed accordingly, and adapted coordinates be 
chosen. The co-frame assumes the structure 



{<) 



( ^ 


N^ 





V 





(11) 



where we let / run horizontally, and /i vertically. The index i now runs from one to three, 
and a is a tangent space index for E4. Analogously, w can be decomposed into S0(3) 
connections T\ :— e^^KL^a^, Ka '■— ^^a ^^^"^ ^^^ rest, i.e., the components of wo- The 
action can now be expressed in terms of the decomposed fields, and it organizes precisely 
as in the case of the electromagnetic fields ^5j : 



5-A 



dt I E-Al 

St 



UG, 



N^C'^ + NC) 



(12) 



= :'tca 



with 



A^^^K + PK, Ef^^detqel, q,, = 4e^,5.,. (13) 

q denotes the metric on St and the dot is the time derivative f^da- We see that A and E 
are conjugate canonical variables. 



{j^,{xlE]{y)] ^ -P5l5]5{x,y). 



(14) 



^ P does acquire physical significance in the case that spinor matter is coupled to gravity, see [14] for details. 



where k := IGttG with G being Newton's constant. The Haniihonian density /lean is a 
sum of constraints, and the constraints 



G. 


= D^^^Ef, C; = E^F: 


, + (• 


■)LG„ 


EfE] 
detE 


'e^hF^,^2{l + /3^)Kl^ 


+ (• 


■)aG.. 



,-.-. , „ (15) 

^' = t; 

By subtracting multiples of Gi from G'^ and C", we obtain a set of independent constraints 
Gi,Ga,C. The constraint equations Gi = Ga — G — 0, together with the evolution 
equations 

{A{x), /iean} = i(.T), {-B(x) , /lean} = E{x) (16) 

are completely equivalent to Einsteins equations. But as it is the case for all reparametriza- 
tion invariant systems, time evolution is a kind of gauge transformation. Concretely, 



{A„G(A)} = -i?(^)A= A 
{E^,G{A)}^[A,E^]= ^ 



g,Ag^ ^ +g,dg^ \ 

geE'^gr' 



e=0 

with (7e = exp(eA) and G(A) = / G^A', so G generates gauge transformations. Moreover, 

{A,CiN)}^C^A, {E,C{N)}^C^E, (17) 

so G generates spatial diffeomorphisms. Lastly G is related to Cnu'^, i.e., to the dif- 
feomorphisms in a direction normal to S. The constraints form an algebra, the Dirac 
algebra, 

{G(A), G(A')} = G([A', A]) {G(A), C{N)} = -GiC^A) {G(iV), C{N')} = G([7V, N']). 

(18) 
The Hamiltonian constraint G is gauge invariant and transforms under diffeomorphisms 
in the expected way, 

{G(iV),G(A)}=0, {G{N),CiN)} = GiCj^N). (19) 

Up to here, the structure is that of an infinite dimensional Lie algebra. But the bracket 
of two Hamiltonian constraints is more complicated, 

{C(N), C(M)} = -'^C(S), S" = ^^(NdbM - MdbN). (20) 

4 dctg 

It contains a function of the phase space point on the right hand side, hence the structure 
is not that of a Lie algebra anymore. 

Before we finish this section, some remarks are in order: 

1. We have just described a canonical formulation of Einstein gravity in Z? + 1 = 4 
dimensions in terms of a phase spaces that is identical to that of SO(D) Yang-Mills 
theory. This formalism relies on a coincidence that only happens for D — 3: An 
SO{D) connection has D{D — l)/2 components, whereas a spatial frame has D. If 
they are to be canonically conjugate variables, they have to have the same number 
of components, which restricts to D = 3. But this does not mean that there are no 
other ways to formulate canonical GR in gauge theory variables. For example, it 
has long been known that GR in D+l=3 has such a formulation in terms of S0(2,l) 
connections and their conjugate momenta, see for example |15j for an introduction. 
More recently it has been shown that gravity in Z? + 1 dimensions can indeed be 
formulated in terms of SO{D + 1) Yang-Mills theory variables [T^, but there are 
additional constraints that have to be implemented in order to bring the number of 
degrees of freedom in line with that of GR. 



2. In the following, we will go over from a formulation in terms of S0(3) to one in terms 
of its covering group SU(2). This must be done in order to couple fermions to gravity, 
but it does not change classical or quantum theory much. The biggest change is that 
also representations with half-integer spin will be allowed in the quantum theory. 

3. As remarked before, in the presence of fermions, physical predictions do become 
dependent on the value of /? [TT . Solving dS/6u! = and re-inserting into the action 
gives an effective 4-fermion interaction with /3-dependent strength. The effect is 
however suppressed by the gravitational coupling constant, and is thus extremely 
small. 

4. We have seen that the constraints commute on the constraint surface G — C = 
C — 0. This means that they form a first class system, and thus the constraints 
can be imposed in the quantum theory as operator equations. This is the result of 
imposing the time gauge. Without it, the situation is more complicated, but also 
very interesting [lIlIllllIllEOlEIlEaEllEllES]. 

2.2 Dynamics in the classical theory 

We saw in the last section that we need to choose the classical variables in which we 
would like to formulate our classical theory with care. Whatever classical choice we make 
will also have an influence on the quantum theory that we obtain. In this section we will 
briefly discuss some aspects of the classical dynamics of general relativity and see that 
also here one needs to think carefully how to formulate the classical dynamics because 
each choice will enter crucially in the properties of the quantum dynamics that we obtain 
later on and discuss in section [3l 

The dynamics of classical general relativity is encoded in the canonical Hamiltonian 
introduced in the last section. As we saw the Hamiltonian consists of a linear combi- 
nation of constraints only and the smeared version of the constraints satisfy the Dirac 
algebra shown in equations (fT8|) . (fT9|) and (l20l) . The constraint C generates spatial diffeo- 
morphisms within the spatial hypersurface, the Gauss constraint generates SU(2) gauge 
transformations and the Hamiltonian constraint generates diffeomorphisms orthogonal to 
the hypersurfaces. Note that the latter is only true on shell, that is when the equation 
of motion are satisfied. From the explicit form of the constraints as well as from the 
Dirac algebra we see that the most complicated among them is the Hamiltonian con- 
straint. In particular it is only the algebra of the Hamiltonian constraints that involves 
structure functions instead of structure constants and thus their algebra is no true Lie 
algebra anymore. As a consequence the corresponding quantum operators need to satisfy 
a complicated quantum algebra in order to be implemented without anomalies. A fur- 
ther complication is that in the representation used in LQG the infinitesimal operators 
for the spatial diffeomorphisms cannot be implemented as operators for the reason that 
their finite counterparts are not weakly continuous (see also section [3.3.3l for more details). 
Hence, strictly speaking one is not even able to check whether the algebra of the quantum 
Hamiltonian constraints is correctly implemented. However, in [26^ an anomaly-free quan- 
tization of the Hamiltonian constraint has been introduced and its anomaly freeness has 
been achieved by quantizing the operator in such a way that the commutator annihilates 
spatially diffeomorphism invariant states. This is what we would expect from an operator 
version of C if it exists. In order to obtain the physical sector of the quantum theory 
later on using the Dirac quantization procedure we look for solutions of the constraints 
in the quantum theory. Also here the complicated algebra of the Hamiltonian constraint 
operators makes our life more difficult and although general solutions to the Gauss and 
diffeomorphism constraint can be constructed one has not been successful in the case of 
the Hamiltonian constraint. However, this is not too surprising since otherwise we would 
be able to write down the general solution of quantum gravity and this is not even possible 
for classical general relativity. During the last years a couple of new proposals have been 



introduced for describing the (canonical) dynamics of LQG. Their aim is to reformulate 
the dynamics in a technically different but physically equivalent way in order to extract 
the physical sector out of LQG. One of the first proposals in this direction is the so called 
master constraint program [27] . The idea of the master constraint is to reformulate the 
dynamics of the Hamiltonian constraints in such a way that the algebra simplifies. For 
this purpose one introduces at the classical level a master constraint involving the squares 
of the Hamiltonian constraints 

Classically due to the square of the Hamiltonian M = is equivalent to C{x) = for all x G 
S. The 1/y/q has been introduced because then M is spatially diffeomorphism invariant 
and thus Poisson commutes with C and also because then M has density weight one which 
will be useful for the quantization of the master constraint. The master constraint satisfies 
the trivial algebra {M, M} = and this carries also over to quantum theory where the 
corresponding operator is required to satisfy this trivial algebra as well. Therefore as far as 
the algebra is concerned the Master constraint simplifies the situation. However, classically 
as well as in the quantum theory at the end we are not only interested in solutions to the 
constraints but we also would like to have so called weak Dirac observables. These are 
quantities that are gauge invariant, meaning that they commute with all constraints of 
the theory on the constraint surface. However, on the constraint surface for the master 
constraints we have {/, M} — for all phase space functions /. Therefore it seems that 
using the master constraint we loose the ability to detect weak Dirac observables with 
respect to the Hamiltonian constraint. However, this is not the case because instead 
of using the condition that weak Dirac observables need to commute weakly with the 
constraint we can also require for a weak observable O 

{0,{0, M}}m=o-0 (22) 

This condition is no longer linear in O as before and hence we get a system of non linear 
partial differential equations for the observables. Formally, it seems that the price we 
have to pay for a simpler constraint algebra is a more complicated equation for weak 
Dirac observables. However, whether the condition in (|22|) is indeed harder to satisfy 
needs to be checked in applications of the master constraint program. As far as the 
quantum theory is concerned the master constraint program has advantages with respect 
to using the Hamiltonian constraint as discussed in |27j . The master constraint program 
has been tested in a variety models in a series of papers [28] . As introduced in [27j we can 
also define a so called extended master constraints that involves also the diffeomorphism 
and Gauss constraint 

M ^ / ,.JCix^l±jf^l±f^G^ (23) 

In this case the constraint algebra trivializes completely and condition ([^^ selects weak 
Dirac observables for all constraints at once. The extended master constraint has for in- 
stance been used in |29j where the semiclassical properties of the quantum dynamics have 
been analyzed in the context of Algebraic Quantum Gravity (AQG). 

Apart from solutions to the constraints and observables at the end of the day we also 
would like to describe the evolution of those observables in order to be able to describe 
dynamics in the physical sector of the theory. Certainly, their evolution cannot be de- 
scribed by the canonical Hamiltonian for the reason that observables, by definition, need 
to Poisson commute with all constraints. That evolution for observables is frozen when 
one considers the canonical Hamiltonian as the generator of their dynamics is known as 



the problem of time in general relativity. Therefore the question arises whether we can 
reformulate the classical dynamics of general relativity in a way that we can improve our 
situation twofold: 

• we would like to get a less complicated algebra for the constraints and 

• we would like to work in a formalism where the evolution of Dirac observables is 
naturally implemented 

The second point can be addressed in the context of the relational formalism introduced 
in [30] and mathematically improved in |31] . The idea of this framework is to introduce 
reference fields that are used to construct Dirac observables and with respect to which the 
dynamics of the remaining degrees of freedom can be described. In the case of GR this 
means to introduce reference fields that label points in space and time when we want to 
construct observables with respect to the Hamiltonian and diffeomorphism constraint. It 
turns out that by a suitable choice of reference fields we can also make progress regarding 
the first point and simplify the constraint algebra. In the following of this section we 
will discuss the classical theory of two particular models where this philosophy has been 
used and the quantization has been performed using LQG techniques. The quantization 
of these models will be discussed in section 14.21 and 14.31 respectively. 

Before going into the details of those models we will briefly comment on two different 
strategies to quantize systems with constraints. Given a classical system with constraints 
we have the option to solve the constraints classically and work with the so called reduced 
(or physical) phase space or we can follow a procedure introduced by Dirac and quantize 
the kinematical phase space and then implement the classical constraints as operators. 
Following the latter we obtain as an intermediate step a kinematical Hilbert space 'Hkin 
and require for physical states that they are annihilated by all constraint operators, which 
are the quantum Einstein equations in the case of general relativity. Solutions to the quan- 
tum Einstein Equations live in the physical Hilbert space Uphys ■ In the case of reduced 
phase space quantization one needs to quantize not the kinematical algebra but the algebra 
of observables and obtains directly the physical Hilbert space Hphys- The dynamics, as 
we will see in the following two sections, is described by a so called physical Hamiltonian 
at the classical level, whose explicit form depends on the chosen model. The Hamiltonian 
is called physical here for the reason that in contrast to the canonical Hamiltonian it is 
not vanishing on the constraint surface. In the reduced case the quantum dynamics is 
encoded in the operator corresponding to the physical Hamiltonian and the quantum Ein- 
stein equations are given by the quantum evolution of the quantized observables generated 
by the physical Hamiltonian operator. More details and concrete applications will follow 
in section 14.21 and section 14.31 Whether one chooses reduced phase space or Dirac quan- 
tization is rather model dependent. In general relativity the constraints are complicated 
and general solutions to the constraints are difficult to construct. That is why often Dirac 
quantization is favored. Furthermore, in general the algebra of observables has a more 
complicated structure than the kinematical one and hence to find representations of it can 
become a challenging task. On the other hand finding solutions that are annihilated by 
all constraints operators and construct the inner product for the physical Hilbert space 
can be very difficult as well and without the physical Hilbert space the physical relevance 
of the results is rather hard to evaluate. The models discussed in the next sections have 
both the properties that the algebra of the observables is isomorphic to the corresponding 
kinematical algebra and thus representations can be easily found. They also share the 
property that not all constraints are reduced at the classical level but part of them are 
solved in the quantum theory. Therefore in both models a combination of Dirac and re- 
duced phase space quantization is used which is also a choice one can make for quantizing 
systems with constraints. 



2.2.1 Brown-Kuchar model 

The idea of the Brown-Kuchar Model [31] is to introduce additional matter dust fields 
that can serve as a reference system for general relativity. In the language of the relational 
formalism we need one reference field for each constraint that we would like to reduce at 
the classical level. The action introduced in [32] has the following form 

S^ust ^-\ I d^X^\Aet{g)\p{g^''U^U, + 1) (24) 

Here M denotes the spacetime manifold, g the spacetime metric, p is the dust energy 
density and C/^ = — T/i + WjS^^^ denotes the dust four velocity, that is itself expressed in 
terms of 7 scalar fields T, S^ , Wj where j runs from 1 to 3. The action is understood as 
a functional of the metric g and the eight scalar field p, T, S^ , Wj that are considered in 
addition to the gravitational and possible other standard model matter degrees of freedom. 
Here we will restrict our discussion to gravity and dust only, but the generalisation to 
additional standard model matter is straightforward. That the action in ([24]) is associated 
with dust is justified by looking at its physical interpretation. The energy momentum 
tensor T^j, is that of pressureless dust. The Euler-Lagrange equations for the scalar fields 
yield that C/^ = g^^U^, is a geodesic congruence and the fields Wj,S^ are constant along 
those geodesies. Furthermore T defines proper time along each geodesic. The canonical 
analysis of the action above shows that the system including the dust involves second 
class constraints. Introducing the associated Dirac bracket and solving those second class 
constraints strongly yields that neither p nor Wj are independent variables but can be 
expressed in terms of the remaining phase space variables. Moreover, as long as only 
functions of T^S^,A and their conjugate momenta P,Pj,E are considered, the Dirac 
bracket reduces to the Poisson bracket, and it is those functions that we are interested in 
the following. For the reason that the fields S^ are constant and the field T defines proper 
time along each geodesies they are a natural choice as the reference fields for space and 
time respectively. The field T will be used to construct observables with respect to the 
Hamiltonian constraint C and the three fields S^ will serve as reference fields for the spatial 
diffeomorphism constraints Ca ■ The choice of reference fields in the relational formalism 
is completely arbitrary and those fields do not even need to be additional matter fields. 
For instance one could also choose gravitational degrees of freedom. However, what will 
be influenced by that choice is the form of the algebra of the constructed observables as 
well as the form of the physical Hamiltonian. Since we are interested in quantizing the 
reduced theory later on we would like to choose those reference fields that lead to a simple 
observable algebra as well as to a physical Hamiltonian, which can be implemented as an 
operator later on. As we will discuss now the dust fields T, S^ satisfy both requirements 
and can be interpreted as a free falling observer that is dynamically coupled to the system. 
The particular choice of T and 5*^ should not be understood as the only convenient clock 
and rods for general relativity but a possible choice that has been made in this model 
and that allows, as we will see in section 14. 2[ to complete the quantization program. 
Likewise to the situation when one chooses a particular gauge fixing it might be necessary 
to consider more than one choice of reference fields in order to consider dynamics in a 
large region of spacetime. By considering the phase space whose coordinates are given by 
T, S^ , A and their conjugate momenta we work in an extended phase space picture where 
time and space are themselves configuration variables. The (first class) constraints of the 
system gravity plus dust have the following form: 



C*°* = C + C"^"''', with C'^""* = -y'p2 + qab(pT^^ + PjSJa)(PT,b + PjSJb) 
C*°* = Ca + Cf''* with C^''''' = PTa + PjS% (25) 

where C, Ca denote the gravitational parts of the Hamiltonian and diffeomorphism con- 
straints shown in (1151) and the Gauss constraint is unaltered the one in (fTSl) . What Brown 



and Kuchar observed in their seminal article was that (i) the fields P,T enter into C'°* 
only in the combination that also occurs in C^"*'*, and using C*°' = those terms can be 
replaced by —Ca and (ii) the constraint C'°* can be solved for the dust momentum P and 
C^°* for Pj . As a consequence one can write down an equivalent set of constraints given 

by 



C'tot ^ p ^ /j ^^^Yl h{A, E) := y/C^ - q<'''CaCb 

Cf = Pj+hj with hj{T,S\A,E) = S^{Ca-hTa) (26) 

where g"'' is understood as a function of E and we had to assume that the inverse of 
S^^ denoted by 3°" exists. In regions of the phase space where det{S^g^) = those fields 
would not provide a good choice of spatial reference fields. In the case of the Hamiltonian 
constraint one obtains a quadratic equation in P and chooses that solution of the possible 
two that yields a positive physical Hamiltonian, which has the correct flat spacetime limit. 
Note that the function h in contrast to hj in (P5|) does not depend on the dust degrees of 
freedom anymore. If this happens one calls a constraint deparametrized, and given this, 
the construction of observables technically simplifies as well as deparametrization for the 
Hamiltonian constraint ensures that the final physical Hamiltonian will be time indepen- 
dent. The constraints in (|26p mutually commute up to SU(2) invariant combinations of 
Gauss constraints 

{&°'{x),C'°\y)}G,^o = {Cf*(x),Cr(y)}G,=o = {C*°*(x),Cf*(y)}G,=o = (27) 

which can be seen by direct computation or by the abstract argument that those con- 
straints are linear in the dust momenta. Since their first class property has not changed 
by writing them in equivalent form, the Poisson bracket of any two constraints can only be 
a linear combination of those. But since they are linear in the momenta and the derivative 
in the Poisson brackets cancels the momentum, the only possible coefficients that are al- 
lowed on the righthand side are zero for C*°' and C,'°* and non-zero coefficients for SU(2) 
invariant combinations of Gauss constraintqj. In contrast the constraints in (P5|) satisfy 
the Dirac algebra and thus by solving for the dust momenta we obtain a simplification 
of the constraint algebra. The construction of observables is performed in two steps now: 
First we reduce with respect to the spatial diffeomorphism constraint by using the fields 
S^ as rods and obtain for the remaining degrees of freedom the following expressions [33] : 

A^{a\t) := f d^x\detiSi)\d{S''{x),a'^)Ais;;ix) 



&j{a',t) := / d^x\dct{Si)\6{SHx),<j')E-jSi{x) 
f(cr^i) :== [ d^x\det{S%)\SiS\x),a'')T{x) (28) 



where we introduced a capital J for the su{2) index in order to distinguish between the 
Lie algebra index and the j — index associated with the dust fields. The interpretation of 
the observables in (P5)) is that they give the value of the fields A,E,T when the dust fields 
S^ take the values a^ . The abstract points x in the spatial manifold S over which the 
integration is performed above have therefore been labeled by the dust fields S^ . Analogous 
observables for the ADM variables where constructed in [32] , however the observabels with 
respect to C*°* where not constructed in [55] but a formal Dirac quantization was used. 
This quantization is formal in the sense that similar to the Wheeler-DeWitt equation 
no representation for general relativity formulated in terms of ADM variables has been 
found so far in which the constraints can be promoted to well defined operators. For the 



When working with ADM instead of Ashtekar variables the Gauss constraint is absent and in this case the 
constraints in (|26|) mutually commute strongly. 
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construction of observables with respect to the Hamihonian and difFeomorphism constraint 
we take the expression for AHa'^jt) and E-'j{a'',t) above and insert them into the standard 
formula for observables, that is a power series in powers of the clock field, in our case T 
and multiple Poisson brackets between the Hamiltonian constraint and AH(T'',t) and 

E-'j {(J^ , t) respectively. Note that only the function h in C'°* will contribute in those 
Poisson brackets, because P commutes with A'Ua^ ,t) and Ej{(t'', t). Explicitly we obtain 



OO -. 

A/(a^r) = exp(x.J = ^-{/i(r),i/(a'=)}(„) 

ES(a^r) = eMXh.) = f2^{Hr),E^j{a'^)}in) (29) 



n=Q 



with 



h{T) := / d^cr(r - T{a))h{a) (30) 

5 

here Xh^ denotes the Hamiltonian vector field oi Ht , S denotes the range of a and is 
also called the dust space and the iterative Poisson bracket is defined as {/, g}(o) ~ 9, 
{f,g}(n) '■= {/j {/)fl'}(n-i)}}- The interpretation of the quantities in (l29l) is that they give 
the values of A and E respectively when the clock field T takes the values r and the rod 
fields S^ take the values a^ . t and a^ can be understood as the physical time and space 
parameter. One can check by explicit computation that A„^((t'^,t) and Ej(cr'^,T) indeed 

Poisson commute with C*°* and C*°' and are thus Dirac observables with respect to the 
Hamiltonian and spatial diffeomorphism constraint. 

The remaining Gauss constraint will be solved in the quantum theory via Dirac quan- 
tization. Once the observables in ((29)) are constructed we are interested in the generator of 



their dynamics. By construction they Poisson commute with C*°' and C*°' and thus the 
canonical Hamiltonian density /lean in ([121) . that generates only gauge transformations, 
cannot be the generator as otherwise the dynamics of the observables would be frozen. It 
turns out that the Hamiltonian density of the generator of the physical dynamics is given 
by the observables associated with the function h in (I^B)) . The observable of a function 
of A, E can be easily constructed due to an automorphism property that the observable 
maps satisfy and hence we obtain 



C2(a,T)-qJ-^-C,Cfc(a,r) (31) 

where C((T, r) and Cj(tT, r) are the observables associated with the gravitational part of 
the Hamiltonian and spatial diffeomorphism constraint that can be obtained by replacing 
A,E by their corresponding observables A, E. Note that although Cfe((T, r) and C(i7, r) 
depend on r the final expression for h does not because on the Gauss constraint surface 
we have 

^ {h(a),h(a')}=0 (32) 

This follows from the fact that C'°*(x) Poisson commutes with itself and since h does 
not include any dust degrees of freedom it immediately follows that {/i(x), /i(2/)}g' =o = 0. 
The physical Hamiltonian is then given by 

Hphys = j dVh(a) = j <fa^lc^~^f^C~Ck{a) (33) 

5 5 

where as before q*-' is understood as a function of E. The equation of motion for the 
observables have the following form 

^ = {A/,HphyJ, ^={E^„Hp,yJ (34) 
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and can be understood as a, with respect to the Hamiltonian and spatial diffeomorphism 
constraint, gauge invariant version of Einstein's equations for the reason that all terms in 
(|34l) are manifestly gauge invariant. In contrast to i/can in the ordinary Einstein equations 
Hphys is nonvanishing on the constraint surface because only the total constraints including 
the dust degrees of freedom are vanishing there. Furthermore Hphys, being itself an 
observable, has the following symmetries on the Gauss constraint surface 

{Hphys, Cj(cr, r)}Gj=o = and {Hphys, h(cr)}Gj=o = (35) 

Those symmetries of Hphys will become important when the quantization of the Brown- 
Kuchar model in section BT^ is discussed. 



2.2.2 Scalar field model 

In the last section we discussed the classical Brown-Kuchar model and saw that the four 
dust fields can be used as reference fields to construct observabels with respect to the 
Hamiltonian and diffeomorphism constraint. The scalar field model, originally introduced 
by Rovelli and Smolin 35 and studied further by Kuchar and Romano [33], considers 
gravity coupled to a massless scalar field whose action is of the form 



S^- 2 



^Jd^Xy^g^^'',p,^cP,, (36) 

M 



The reason that one considers a massless scalar field is that otherwise the resulting Hamil- 
tonian constraint would not be deparametrizable having the consequence that we end up 
with a time dependent physical Hamiltonian. Here we introduce only one additional mat- 
ter field coupled to gravitjO and thus we are only able to reduce one of the constraints 
at the classical level. That will be the Hamiltonian constraint. The remaining diffeomor- 
phism and Gauss constraints will be solved in the quantum theory in this model. One of 
the motivations for this model is that in recent models in loop quantum cosmology (LQC) 
also a scalar field is chosen ^7\ as a clock, so this model here could be understood as 
the generalization of the APS model in [37] to full LQG. Further motivations will become 
clear in section IT2] and l473l when the quantization of the scalar field and the Brown-Kuchar 
model is discussed. A 3-t-l-split of the action above yields the following total first class 
constraints 




v/det(q)q^V,a0. 



, v/dct(q) 

C = Ca + Ct with C^^^' ^ TT(f>^a (37) 

where as before C, Ca denotes the gravitational parts of the Hamiltonian and diffeomor- 
phism constraints shown in (|15p and n denotes the momentum conjugate to (f>. As before, 
the Gauss constraints is not affected by the scalar field and thus it is still the one shown in 
(J15p . The Brown-Kuchar mechanism, that is replacing the terms that involve derivatives 
of (/) by C^ and this then again by —Ca and then solving the Hamiltonian constraint for 
the scalar field momentum n can be applied. One obtains a fourth order polynomial in 
n and out of the four solutions only two are non trivial in the cosmological sector. The 
remaining ± sign ambiguity defines different regions in phase space and we choose tt = +h 
for the reason that only this solution includes homogenous and isotropic models. Using 
this choice one obtains an equivalent Hamiltonian constraint of the form 



C*°* = TT-h with h{A,E)■.= ^J-y/^C + y/^^/C^-q<'>'CaCb (38) 



Note that also in the scalar field model the generalization to gravity plus arbitrary standard model matter 
is straightforward. 
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Now we are in an analogous situation to the Brown-Kuchar model in section [2.2.11 We 
realized that also the system of gravity coupled to a massless scalar field deparamctrizes 
and we could proceed as before. This would mean to use the scalar field cf) to construct 
Dirac observables with respect to the Hamiltonian constraint. What we would obtain are 
quantities A^ (r, x) and £'" (r, x) that commute with C'°* in (|55|) but not necessarily with 
the diffeomorphism or Gauss constraint. The latter two will be solved in the quantum 
theory via Dirac quantization. The algebra of A^(t, x) and E'^{t, x) is isomorphic to the 
one of kinematical quantities and their dynamics, that is evolution with respect to the 
physical time parameter r, is generated by a physical Hamiltonian of the form 

H^i^ys^ I d?xh{A,E){x) (39) 

E 

with h{A,E) given in equation ([38)) . Alternatively, we cannot construct those observables 
with respect to the Hamiltonian constraint at the classical level but in the quantum theory 
and work with so called quantum Dirac observables, whose classical limit corresponds to 
those observables mentioned above. This will be explained more in detail in section 
14.31 In this case all constraints will be solved by Dirac quantization in the quantum 
theory. Likewise to the Brown-Kuchar model, also here the Hamiltonian densities strongly 
commute and iJphys is invariant under spatial diffeomorphisms, that is 

{Hphys, h{x)}Gj=o = and {i^phys, Cair, x)}g,=o = (40) 

In contrast to the Brown-Kuchar model a quantization of this partially reduced system will 
not yield the physical Hilbert space. The quantum theory of this model will be discussed 
in detail in section 14.31 Finally, let us mention that there exist another model introduced 
in [38j where only one instead of four dust fields are considered. This model seems to be 
a special case of the Brown-Kuchar model, where the momentum density of the dust has 
chosen to be equal to zero. 

3 Quantum kinematics 

This section includes a introductory discussion on the quantum kinematics of loop quan- 
tum gravity corresponding to the sector of the theory where the constraints have not been 
solved yet. For the reason that the kinematical representation is rather different from the 
usual Fock representation we introduce it in section 13.11 and discuss its properties. Af- 
terwards in section [3?2l we discuss geometric operators in the framework of loop quantum 
gravity. These are operators that correspond to classical geometrical objects like length, 
area and volume. Section 13.31 deals with the question how the classical Gauss and dif- 
feomorphism constraints can be quantized and how solutions to the constraint operators 
can be constructed. Finally, in section 13.41 we present more details on the kinematical 
representation and explain how it can be generalized. 

3.1 Ashtekar-Lewandowski representation 

We will now quantize gravity according to the algorithm for the quantization of con- 
strained systems devised by Dirac (for the original account, see his Lectures on Quantum 
Mechanics, for a modern treatment, see [2]). This means we proceed in two steps. 

1. Quantization of the canonical variables ("kinematic quantization") 

2. Impose the constraints as operator equations on states, and solve these equations to 
obtain physical states. 

The first step is what we will discuss in the present section. What we want is a represen- 
tation of the canonical commutation relations 

[Al{x),E^{y)]^^h(3Sl6]S{x,y) (41) 
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on a Hilbert space. Note that kH = £p, the Planck area. Fields evahiated at points 
are usually too singular to give good operators in the quantum theory. Thus one has to 
form suitably integrated (" smeared" ) quantities that correspond to well defined operators 
in the quantum theory. Poisson brackets then suggest commutation relations for these 
quantities, and one obtains an abstract algebra of operators. We will soon see that details 
can matter in this context. Different choices of smearing can lead to different algebras and 
hence to different quantum theories. In fact, in LQG we make a different choice of algebra 
than is customary in Yang-Mills theory [311 HOI Ell SSI SSI Hi]- In the latter case, both 
the algebra and its representations used in the quantum theory make use of the metric 
as a classical background field in their construction. In general relativity, the metric is 
dynamic and hence can not be used as a background field. Moreover, a splitting of metric 
into background and dynamical part, while very useful in practical applications, is not very 
natural from a fundamental perspective. Hence, the algebra and its representation used 
in LQG does not make use of any background metric. This makes LQG a very unusual 
quantum field theory. To illustrate this, we first consider the case of electromagnetism. 
The usual quantum theory is obtained by declaring 

f rAay/d^d^x, f flE'd^x =ih f rf^^^d^d^x id, (42) 

and by defining the ground state by 

aim - 0. (43) 

Here /, /' are arbitrary smearing functions, and the definition of the annihilation operators 
a makes use of the metric q in various ways. But one could also define 

E{S):^ f E'^eabcdx'Adx' (44) 

Js 

where S is an oriented surface and eabc is the tensor density that is equal to the totally 
anti-symetric symbol in any coordinate system. We note that E"" has density weight +1 
whereas eabc has weight —1, so the integrand is a two- form and the integral, using the 
orientation of S, is hence coordinate independent. Similarly, defining 



Aie)= Aadx'', (45) 

where e is a curve, one finds (by a limiting procedure from the Poisson brackets of the 
point fields), 

[Aie),E{S)]=ihI{e,S) id, (46) 

where /(e, S) is the signed intersection number for S and e. The metric has thus dropped 
out of all definitions and relations. 

A similar thing can be done for gravity. For a surface S and a Lie algebra valued 
smearing field n on S we define 

EniS) := / n'E^eabc dx'' A dx'. (47) 

Js 

For A we use the quantity analogous to exp(iA(e)). We chose a local trivialization and 
define the holonomy 

he :=Pexp / A (48) 

J e 

= 1+1 A(e(t))e"(t)dt+ / dh I dt2Aa{h)e-'{h)Aa{t2)e''{t2) + ... (49) 



Jo Jt 
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which is an element of SU(2), and gives the parallel transport map from the fiber over 
the beginning point 6(e) of the edge to the fiber over of its final /(e), for the chosen 
trivialization. Under gauge transformations, i.e., changes of trivialization, 5 : M — >■ SU(2) 



K^g{b{e))Kg(f{e)r\ (50) 



One finds 

[En{S),h. 



J '-Y-heiTjU^ {p)he2 for a single transversal intersection S* n e — {p} 
~\o ifS'ne = 

(51) 
where {r,} is a basis of su(2), and ip = hn is (a multiple of) the Plank length. It is 
convenient to slightly generalize these variables. Given a graph of paths 7 — {ei, e2, . . . e„} 
and a function / :SU(2)" -^ C, one obtains the functional 

f,[A]:^f{K,[A],h,,[A],...h,,AA]) (52) 

A functional f is called cylindrical with respect to 7 (written / G CyL) if it is of the above 
form, and simply cylindrical if it is of the above form for some graph 7. We note: 

• A given cylindrical functional is cylindrical on many graphs. Consider the example 
of a function f-y[A\ — /(/ie[A]), which is cylindrical w.r.t. the graph 7 — {e}. Now 
consider a second graph 7' = ei,e2,e3, with ei o e2 = e, and ea independent of e. 
Then /^ is also cylindrical w.r.t. 7', as it can be written purely in terms of holonomies 
along edges in 7', f^[A] = /(/ieJ^]/ieJA]). 

• For two cylindrical functions which are cylindrical on graphs with smooth edges, one 
can not always find a finite graph such that they are both cylindrical w.r.t. to that 
graph, because they can intersect infinitely many times. But for more regular edges, 
for example analytic or semi-analytic (roughly speaking: piecewise real analytic [IS]) 
ones this can not happen, and so one can always find such a graph. As a consequence, 
such cylindrical functions are closed under addition and multiplication and thus form 
an algebra, denoted Cyl. In the following, we will always assume edges (and also 
surfaces) to be semi-analytic. 

We can use these observations to write the commutator between the canonical variables 
in a relatively concise form. For this, we assume without loss of generality that a surface 
S and a graph 7 intersect only in vertices of 7. The commutator then reads: 



[f^,ErXS)]^XrXS){f,)^^ ^ n\v) 
where V{'^) denotes the set of all vertices of 7 



^ «(5,e,z;)jJ''-'/ 



(/lei,/le2,---)) (53) 



)0 if e intersects S tangentially in v or does not intersect S at all 
1 if e intersects S transversally and is above S (54) 

— 1 if e intersects S transversally and is below S 

and 

-^v e) . , . . { Lt \ . , , f e ingoing at u 1 ,^^- 

4"'"=^ = id «) id (g) . . . «) <^ „l \®\d® ..., when <^ 7 ■ ^ } ■ 55) 

'^ I "fc J I ^ outgoing at w J ^ ' 

R, L denote the right/left invariant vector field on SU(2) associated with a basis r^ of 
su{2). The additional factor of 1/2 in (15^ as compared to ([ST]) is due to our assumption 
that edges must end on the surface. An edge that continues on both sides of the surface as 
in ([5T|) will thus count as two separate edges in (|53p . For general surfaces, the commutator 
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above may not be a cylindrical function, again because edges and surfaces can interact each 
other infinitely often. Thus we also restrict the surfaces to be in a suitable regularity class. 
Then the operation Xn{S) defined above is a derivation on the space Cyl of cylindrical 
functions. We note that the commutator has the Jacobi property, so 

[f,[E^{S),Er,,{S')]] = [X^{S),Xr,,{S')]{f) (56) 

and the commutator on the right hand side is non-vanishing in general. Thus we find that 
the operators corresponding to the spatial geometry do not commute. 

The objects En{S), together with the cylindrical functions Cyl subject to the above 
commutator relations form the kinematic algebra 21. Since it does not make reference to 
classical geometry on E, diffeomorphisms act in a simple way: 

MfM] := f{cf>*A), a^{En{S)) = E^,^{<p{S)) (57) 

are automorphisms of 21. A similar statement holds for gauge transformations. 

To implement the constraints we need a representation of 21, i.e., a mapping of 21 into 
the operators of a Hilbert space. There are many representations of 21, but one of them is 
special and therefore most important for LQG, the Ashtekar-Lewandowski representation 
of 21 [m |42l |43l |46]. To see how it works, note first that an inner product on Cyl can be 
defined by 



{f-/\fj)-= / dfi{gi) .. .d^i{gn)f {91,92,- ■ ■gn)f'igi,92,- ■ -gn)- 

JSU(2)" 



(58) 



The measure d/i used above is the Haar measure on SU(2), and we have assumed without 
loss of generality that the two functions are cylindric w.r.t. the same graph, as discussed 
below ([5^ . Closure with respect to the corresponding norm gives a Hilbert space Ti- 
lt can be shown that this space has a very suggestive structure, H = L^{A,dfiAL), the 
square integrable functions over a space of distributional connections, with respect to a 
certain measure. 

The action of the basic operators in this representation is analogous to that found in 
the Schrodinger representation of quantum mechanics: 

nifMA] = f[AMA], 7riE^{S)MA] ^ iX^{S)^)[A], (59) 

where we have assumed that ^I* is smooth enough for Xn{S) to act. For example, \1/ 
could be a cylindrical function based on a differentiable function on some power of SU(2). 
But the properties of this representation are very different from those of the Schrodinger 
representation of quantum mechanics, and of the representations encountered in standard 
QFT. For example, eigenstates of the fiuxes, i.e., the momentum variables, are normaliz- 
able, as we will see in a moment. Also, there is a precise analogue of this representation 
in the case of a scalar field, and it is unitarily inequivalent to the standard representation 
for a scalar in flat or curved background [48] . 

The representation has several useful properties. It is irreducible and faithful. No 
background geometry was used in the definitions, so it carries a unitary representation of 
spatial diffeomorphisms and gauge transformations. It has an invariant and cycliqj vector 

n. 

The Hilbert space H has a very useful orthonormal basis. To explain, let us first 
consider a general compact Lie group G. Then there are two natural representations of G 
on Hg — i^(G,d/i), the left- and right-regular representations 

(pl(<?)/)(3') - fig'g-'), {pn{9)f){g') = figg')- (60) 



''Cyclic means that {an|a £ 21} is dense in "H. 
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They both decompose into irreducibles, and since the two representations commute, there 
is a common basis of eigenvectors of the Casimir operators. Let tt be an irreducible 
representation of G, then 

y(7r, m) :— span{7rmn(-)l" = Ij 2, . . . dimTr} is left invariant by pL, (61) 

y(7r, n) :— span{7r„i„(-)|m = 1,2,... dimTr} is left invariant by pR. (62) 

The representation V^n, m) induced by pL is tt itself. The one induced by pR, on V{it^ n) 
is its dual, 7f, i.e. ¥((7) = ■K{g~^)'^ . The Peter- Weyl theorem now asserts that each irrep 
arises in the decomposition of the regular representations, and even more, that their matrix 
elements give a basis of "Hg- Pick, for each equivalence class of irreps of G a representative 
TT, then the set of all Vdim TTTTmn for all equivalence classes form an orthonormal basis of 

na. 

Now back to the LQG setting. Let H^ — Cyl . On the one hand, 'H-y is a subspace 
of H, on the other hand it is isomorphic to i^(SU(2)"). Thus, an orthonormal basis for 
"H^ is given by 



^/{2n + l)7^:{2j^^+l)'nkMhe^[A])■nk,l,iheM]) ■ ■ ■ -^ k„iJK^[A]) ■ (63) 

where the ji, J2, • • -jn label irreducible representations of SU(2). We note, however, that 
in general Jiry JL "H-y'. Given 7, there are functions in "H^ that do not depend on the 
holonomy along some edge of 7, or that only depend on the product of two holonomies 
along adjacent edges in 7. Such functions are also elements oiH-y' for some 7' that has less 
edges and vertices. This makes linear dependencies among the elements of different graph 
Hilbert spaces rather generic. Take for example 7 = {e}, 7' — {61,62} with 6 = ei o 62. 
Then 

TTmniheiA]) = ^ 7r„™/ (/lei [^] )7rm'„(/le2 [A]). (64) 

m' 

Therefore we will introduce a family of slightly modified Hilbert spaces H' , which give a 
decomposition of V. into orthogonal subspaces. But first we need to discuss the transfor- 
mation properties of vectors under gauge transformations. 

We start by considering just a single edge e. With respect to gauge transformations g, 
the vectors 7r^„(/ie) transform under the tensor product j ^ j, and can be visualized as 
the edge with representation j sitting at its endpoint and representation j at its starting 
point. When several edges meet at a vertex v, contractions of the matrix indices of the 
representation matrices at that vertex can be done and correspond to vectors in the tensor 
product 

Uv^i (g) JeJ (g) j (g) Te] (65) 

\ e incoming atv j \ e outgoing atv j 

To give an orthogonal basis of this space, one can simply decompose it into irreps, 

H. = Q /, (66) 

I 

where ci counts the multiplicity of the spin /-representation. When we apply this to the 
situation in LQG we obtain the following decomposition. Given a graph 7, 

^7 = 0^,j = 0^,j,r- (67) 

J IS 

Here we have first decomposed into spaces in which the assignment of irreps to edges 
(labeled by j) is fixed, giving essentially the tensor products of the spaces (|65l) . Then we 
have further decomposed according to (j66p . labeling the irreducible subspace chosen at 
the vertices with /. 
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Now we can go back and remedy the problem that the decomposition into T-L^ was 
not an orthogonal one (following [5]). Given again a graph 7, we can call a labeling j of 
edges and / of vertices with irreps admissible if no two-valent vertex has been assigned the 
trivial representation I = 0, and none of the irreps assigned to the edges is trivial. Then 
we set 

jd admissible 

and obtain the desired orthogonal decomposition 



^=0^;- (69) 



3.2 Geometric operators 



One of the special properties of the representation in LQG introduced in the last section 
is that one can define operators corresponding to geometrical objects such as volume area 
and length. Among those the most simple operator is the area operator from the point 
of view of the construction of the operator as well as with regards to the spectrum of the 
operators. The area operator was first introduced by Smolin j49j and then further ana- 
lyzed by Rovelli and Smolin in the loop representation [SD]- Ashtekar and Lewandowski 
discussed the spectrum of the area operator in the connection representation in [46] . In 
this section we want to discuss the implementation of the area operator as well as its 
spectrum in detail. At the end of the section we will briefly comment on the volume and 
length operator. 

The classical area functional associated to a surface S is given by the following expression 

Ar{S)^ f (fu^det{X*q){u) (70) 



u 
The ADM 3 metric is denoted by (7, A" : [/ — ^ S" is an embedding of the surface, where 
U <zM? and X* denotes the pull back of X. The coordinates on the embedded surface 
S are given by the embedding functions A° with a — 1,2,3 and let us denote the two 
coordinates by which the surface is parametrized by ui and U2- Given the embedding we 
can construct two tangent fields on S 

_ dX" _ dX°- 

and also a co-normal vector field tt-q that is determined from the condition 

naXl^ = Q for i = l,2 (72) 

The determinant in the area functional can be expressed as 

det(A*(j) — <7„i„i<7„2M2 ~ 9«iM29«2Mi — \^%i^ .ui^Z-^^ ,u.^ ~ ^%x^ ,U2-^'!^U2-^ .ui) Qabqcd 

. (^^) 
In order to quantize the area functional we need to express it in terms of Ashtekar variables. 

For this purpose we consider the expression det(g)7T,anhg°'' and use that we can express 

the inverse metric as 

Furthermore we see from ([7^ that Ua — eabcX'^^_^X'\^^ yielding 

dei{q)nanbq''^ = det(g)nanb- . '^'""^'^'"'^qceqdf 

9Mltll9M2"2 '?MiM29'U2"l \'") 
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The inverse metric has a simple form in Ashtekar variables given by q""^ — EjE^S^'' / det(£') 
and depends only on the densitizcd triad. From E^ = y/qe'^ we get det(i?) = det(q) yield- 
ing 

det{q)q'''' = E'^Eld^'' (76) 

from which we can conclude using (|75p that 



^Aet{X*q) = ^nan,,E<^El5^J^ (77) 

Note that often one choses the basis Tj :~ — iCTj/2 in su(2) with aj being the Pauli matrices 
for which the Cartan-Killing metric on su(2) rjjk becomes "qjk '■= Tr(ad{Tj)ad{Tk)) = —2Sjk 
and then one uses the Killing metric in the expression above and adjusts the pre-factors 
accordingly. The area functional is completely expressed in terms of Ashtekar variables 
now. In order to quantize the area functional we need to choose a regularization of the 

n 

classical expression. We choose a partition of U into closed set such that U — [J Ui. 

1=1 
Then we can express the area functional as a sum of integrals over the individual sets Ui . 

Furthermore, we assume that the area of each Ui is e^ and we will denote a point centered 

in Ui by v. 

Ar{S) ^ J2 I d?u^E'^El5J'^nanb{X{u)) 



Y,^^\JE';El53kn,Tn,{v) 

i=l 
n 

Y,J{e^E-naWElnb)5i^^{v) 



= J2yEj(Su,)E,{Su,)6iHv) (78) 

i=l 

In the limit n — > cxd the Riemann sum above yields exactly the area functional. As can 
be seen in the last line we managed to express Ar{S) as a function of the classical fluxes 
Ej{Sui) — Er-iSui) for which well defined operators exists. For this reason we obtain 
the quantum area operator simply by replacing the classical fluxes by their corresponding 
operators. Recall that the flux operators can be written in terms of right and left invariant 
vector fields R^ and L^ associated to the edges of a given spin network function. 

e at V 

The term k,{S, e, v) is as before +1 and —1 respectively for edges that intersect S transver- 
sally and are above S and below S respectively at v and zero in all other cases. For the 
area operator we obtain 



-'\ 



' Y, n{Su.,e,v)Ji^.^))\jr (80) 



If the partition is chosen fine enough so that only one intersection point exists in Sjji 
instead of summing over all sets Ui in the partition we can sum over all intersections 
points of an edge of type up or down. Note that in each intersection point more than 
one edge could meet. In this way we will loose all surfaces Sui, which do not contain 
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an intersection point of an edge of type up or down. However, these would anyway not 
contribute since for them K(e, Su) = 0. Thus we can reexpress the area operator as 



MS)T^jj-^ E JiT.<Sui.),e,v)J(^^^))\jj 



where the set P{S) of intersection points of edges up and down is given by 
P{S) ^{veeCi S\k{S, e,v)^0, e G E{j)} 



(81) 



(82) 



Let us now discuss the spectrum of the area operator. At each intersection point x we 
have edges of type up, edges of type down and edges of type in that will not contribute to 
the spectrum. In order to write the expression under the square root in (|8ip in compact 
form we introduce the following operators: 



j: 



eeE{v.u) 



(v,e) 



J. 



v.d 



eeE{v,d) 



{v,e) 



(83) 



Here E{v,u),E{v,d) denotes all edges of type up and down respectively that intersect 
each other in the point v. Then we have for each intersection point v 



( 






E k(S', e, u)J, 



{v.e) 



, eGB(7) 



fjv,^ 



J 



V.d 



2(J"^")2 + 2{J"-'^f - (J"'" + r^'^f (84) 



We used in the second line that [Jl'^jJ^' ] = 0. Furthermore the operators (J"'")^, 
(^jv,dy, ^^^ (^J""^^ + j^'"*)^ mutually commute. Moreover, we choose an explicit basis tj = 
—iaj/2 for which the operators J^^-^' satisfy the usual angular momentum algebra given 
by [if ■"\ JJ"'"^] = e»,fc4''"^ Then we have that the operators (j(^'''))2 = 6^^ jf^^W^^^"^ 
locally act as 



5'^R,R, =-{R\R)= -Asu(2), or - 6'^ L,L, = -{L \L) = -A 



SU{2} 



(85) 



where —Asu(^2) is the positive definite SU(2) Laplacian with spectrum j{j + 1), due to 
our choice of basis for su(2). Hence the same holds for the operators (J^'")^, (J"''*)^, 
and (J"'" + J"''')^, they act as Laplacians in the respective direct sum of representations. 
Therefore the spectrum of the operators involved in (|84l) can be easily computed and we 
obtain 

Spec{Ar{S)) = ^ ^ 
veP{s) 



'2juAJu,v + 1) + 2j„,„(j„,„ + 1) - ju+d.viju+d.v + 1) (86) 



Here ju,v,jd,v denote the total angular momentum of the edges of type up (down respec- 
tively) at the intersection point v and ju+d,v total coupled angular momentum of the up 
and down edges whose values range between \ju,v — jd,v\ < ju+d,v < ju,v + jd,v- Let us 
consider the eigenvalue at one intersection point v. The smallest possible eigenvalue that 
we can get occurs when either ju,v = and jd,v = ^ or vice versa. The eigenvalue denoted 
by Ao is non vanishing and given by 



Ao = 11(3^ 



(87) 
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and is known as the area gap in LQG. The area gap plays an important role in the 
description of black hole physics within LQG and black hole entropy calculations can be 
used the fix the value of the Immirzi parameter /3. 

Finally let us say a few words about the other two geometrical operators, the volume 
and length operator in LQG. The volume operator enters crucially into the construction 
of the dynamics of LQG for the reason that the classical co-triad is expressed as the 
Poisson bracket between the connection and the classical volume functional using the 
Thiemann trick (see section [4.ip . The volume operator can be quantized in a similar 
manner than the area operator because the classical three dimensional volume of a given 
region R is the integral over R of ^det{q) whose expression in Ashtekar variables is just 



det(£''')|. Thus it is again a function of the electric fields only and can be after a 

suitable regularization again be expressed in terms of flux operators. In the literature 
exist two different volume operators one introduced by Rovelli and Smolin (RS) [50^ and 
one introduced by Ashtekar and Lewandowski (AL) [JT], which come out of a priori equally 
justified but different regularization techniques. Both volume operators act non-trivially 
only on vertices where at least three edges intersect. At a given vertex the operators have 
the following form 




Crs 2^ \/\Qijk\ 

einejCiCK—v 



^ e{ei,ej,eK)QijK\ (88) 



eiHejOeK—v 

Here Qijk '■— ^ijkJl^JljJeK ^^ ^^ operator involving only fiux operators and thus right 
and left invariant vector fields and CrsjCal are regularization constants. The sum runs 
over all ordered triples of edges intersecting at the vertex v. The main differences be- 
tween these two operators is that the RS-operator is not sensitive to the orientation of the 
triples of edges and therefore also planar triples of edges will contribute. The AL-operator 
has likewise to the K{S,e,v) in the area operator a similar sign factor e(e/, ej, e^) that 
vanishes whenever the triple of edges Ci, Cj, Ck intersecting at a vertex v are linearly depen- 
dent. Furthermore the sum over triples of edges involved in both operators, occurs outside 
the square root in case of the RS and inside the square root in case of the AL-operator. 
The spectral analysis of the volume operator is more complicated than for the area operator 
and can in general not be computed analytically. A general formula for the computation 
of matrix elements of the AL-volume operator has been derived in |52j . Those techniques 
have been used to analyze the spectrum of the volume operator numerically up to a ver- 
tex valence of 7 in a series of papers , 53| . Their work showed that the spectral properties 
of the volume operator depend on the embedding of the vertex that enters via the sign 
factors e(ei, Cj, e^) into the construction of the AL-operator. Particularly the presence of 
a volume gap, that is a smallest non vanishing eigenvalue, depends on the geometry of the 
vertex. For a more detailed discussion about the spectrum of the volume operator we refer 
the reader to the article of Brunnemann in these proceedings. A consistency check for 
both volume operators has been discussed in [54] where the Thiemann trick has been used 
to define an alternative flux operator. Those alternative flux operator is then compared 
to the usual fiux operator and consistency of both operators could for instance fix the un- 
determined regularization constant Cal — ^p/V^S in the volume operator. Furthermore, 
the RS-operator did not pass this consistency check and the reason that it worked for the 
AL-operator is exactly the presence of those sign factors e(e/, ej, ex) in the AL-operator. 

A length operator for LQG was introduced in [55]. The length operator is in some sense 
the most complicated one among the kinematical geometrical operators. Let us recall the 
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the length of a curve c : [0, 1] ^- S classically is given by 



e{c) = f .Jqat{c{t))c-{t)d'>{t)dt = f ■sJelic{t))ei{c{t))d-it)db{t)S.jdt (89) 



here c" denotes the components of the tangent vector associated to the curve. When we 
express the metric qab in terms of Ashtekar variables we obtain 

Qab - -^<^acd<^befiijk<^ det(E) ^ ^ 

which is a non-polynomial function in terms of the electric fields and therefore a regular- 
ization in terms of flux operators similar to the area and volume operator does not exist. 
Furthermore the denominator being the square of the volume density cannot be defined 
on a dense set in "Hkin because it has a huge kernel. One possibility to quantize the length 
used in ^5S] is to use for the co-triads that occur in (15^1) the Thiemann trick and replace 
them by a Poisson bracket between the connection and the volume functional. This yields 
a length operator that involves a square root of two commutators between holonomy op- 
erators along the curve c and the volume operator. In this way the inverse volume density 
can be avoided and the volume occurs only linearly in the commutator. Also, the length 
operator does not change the graph or the spin labels of the edges likewise to the area 
and volume operator. However, since the length operator becomes even a function of the 
volume operator its spectral analysis becomes even more complicated than for the volume 
operator itself and very little about the spectrum of the length operator is known except 
for low valence vertices. 

Another length operator was introduced in [56| where the Thiemann trick was not used 
for the quantization. The regularization adapted in [5 6) is motivated from the dual picture 
of quantum geometry and uses that the curve can be expressed as an intersection of two 
surfaces. This allows to express the tangent vector of the curve in terms of the normals 
of the surfaces. The inverse volume issue discussed above is circumvented by using a 
Tikhonov regularization for the inverse RS- volume-operator. For this length operator the 
spectral properties have only be analyzed for a vertex of valence 4, which is monochro- 
matic, that is all spins are identical. Recently an alternative length operator for LQG has 
been discussed in [57, where a different regularization has been chosen such that the final 
length operator can be expressed in terms of other geometrical objects the area, volume 
and flux operators. In this work the AL-operator is used and the inverse volume operator 
is also defined using a Tikhonov regularization similar to the one in [56) . 

3.3 Solving Gauss and diffeomorphism constraints 

In accordance with Dirac's program for constraint quantization an implementation of the 
classical constraints as operators on the kinematical Hilbert space is needed. In section 
I3.3.1l we will discuss this quantization for the Gauss constraint and explain how solutions 
to the constraint can be constructed. The quantization and the construction of solutions 
to the spatial diffeomorphism constraint requires a framework called refined algebraic 
quantization that will be introduced in section 13.3.21 and applied to the diffeomorphism 
constraint in section [3.3.31 

3.3.1 Quantization and solution of the Gauss constraint 

The classical expression for the smeared version of the Gauss constraint is given by 

G(A) = I Sxk\DaE^) ^ - f d^x{DaA')Ej (91) 
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This expression involves the densitized triad integrated over a 3-diniensional integral with 
smearing function (DaA.^). When the holonomy flux algebra is computed by means of 
a regularisation of holonomies hl{A) and fluxes E'^ {S) both the holonomies and fluxes 
are regularized by expression that involve three dimensional integrals and e and e' denote 
regulators. For this reason one can follow the same regularization strategy for the compu- 
tation of {hf,{A), G{A)} with the difference that we do not need a regulator for G(A) since 
it involves already a 3 dimensional integral. Computing {hl{A), G(A)} and removing the 
regulator afterwards yields 

1 

{he{A),G{A)} - ^^ f dt{DaA^){e{t))e-{t)he{0,t){A)T,he{t,l){A) (92) 



In the following we will choose an explicit basis for su{2) given by tj := —iaj/2 with 
(Tj being the Pauli matrices. In this case the structure functions fjk£ are given by the 
Levi-Cevita symbol fjM — ^jm- The expression /ie(0,t) satisfies the differential equation 

^^^^ = -/ie(0,i)A(e(t)) with A{e{t)):^A^Mm''{t)r, (93) 

The holonomy is defined as the unique solution he{A) :— /ie(0, 1) with /ie(0, 0) = ls£/(2)- 
The covariant derivative in equation (|92|) is phase space dependent since it involves the 
Ashtekar connection A^^. However, as we will show the expression under the integral can 
be expressed as a total time derivative of holonomies and A := A^Vj. As a first step we 
consider 

(D,A)e»(t) 



dk 

dt 


+ eM^aASe 


'^{t) 


dk 
dt 


+ Alk'[TkM 


e\t) 


dk 
dt 


+ [AA] 





(94) 

Furthermore we can use the product rule for the term ^ (/ie(0, t)A/ie(i, 1)) and for the 
individual terms using the differential equation in (|93p 



to obtain 



-/ie(0,i) = K[Q,t)A{e{t)) 

jh,{t,l) = -A{e{t))K{t,l) (95) 



^ (/le(0, t)kh,{t, 1)) - /ie(0, t) (^ + [A, A] ) h,{t, 1) (96) 



Consequently, we obtain for the Poisson bracket the following expression 

{K{A), G(A)} = ^ (kiemhM) - KiA)kieim (97) 

Knowing how G(A) acts on holonomies we can generalize the action from holonomies to 
cylindrical functions and express it again in terms of right and left invariant vector fields. 

{/,(A),G(A)} ^ ^ J2 (A^"(6(e))i?|-A^(/(e))L^^)/,(A) 






^ A^(«)( Y. <S,e,v)J^''''^f^)iA) (98) 
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where we rearranged the sum over edges as a sum over vertices and sums of the edges 
startmg or ending at each vertex yielding k{S, e, v) = +1 and k{S, e, v) — +1 respectively. 
We realize that action of G(A), likewise to the action of the flux, can be expressed as a 
vector field on C°°{A). The quantization of G'(A) can then be obtained as an extension 
of G(A) to cylindrical functions e C°°{A) where A is the quantum configuration space of 
generalized or also called distributional connections. 



j«^2 



G(A) = ^n[/^(A),G(A)] = ^ ^ A^(«)( J] «:(^,e,z;)jf ^^Vt) (^) 

MA) (99) 



^ E A^(-) T. R^- E L) 



2 

wev(7) 



I eeB(7) eG-E(7) / 

\v=b(e) v=f[e) I 



where we used in the last line the explicit expression for J and K{S,e,v). Solutions 

to the Gauss constraint G{A)fj(A) — are those cylindrical functions for which the 
ingoing and outdoing edges at each vertex couple to the same resulting angular momentum 
Jin = Jout so that the total angular momentum J = Jin €5 Jout =0. In terms of the 
decomposition of the kinematical Hilbert space in terms of interwiners the gauge invariant 
Hilbert space is just the subspace where one associates to each each vertex an intertwiner 
that projects on the trivial representation. That this constructs gauge invariant spin 
networks functions follows from the fact that the holonomy transforms under a finite gauge 
transformation as he{As) = g{b{e))he{A)g~^ {f (e)) , where g G SU{2) and 6(e) and /(e) 
denote the beginning and final point of the edge, and hence the gauge transformations 
act at the vertices of the spin network function only. Furthermore, at the vertices the 
spin network functions transform in the resulting representation that the individual edges 
couple to, thus 

^Sn = 0^7j,r=o (100) 



3.3.2 Refined algebraic quantization 

In this section we briefly review reflned algebraic quantization that provides a framework 
to solve first class constraints in the quantum theory. The application of reflned algebraic 
quantization in the context of Dirac quantization has been analyzed by Giulini and Marolf 
in |58[ I59j . A review article about reflned algebraic quantization and group averaging can 
be found in [50] . 

In order to formulate the flnal quantum theory we are looking for solutions to the con- 
straints in the quantum theory, hence we want to flnd states that are annihilated by the 
constraint operators, that is 'Gf/'phys = 0'. In general the constraint operators G are 
complicated functions of holonomies and fluxes and we are looking for the value zero in 
their spectra. In the simplest case the operators that we need to consider are bounded 
and their spectra is discrete. Then the operators are deflned on the whole Hilbert space 
"Hkin and the eigenvectors of the operators are elements of 'Hkin- As we will see later on 
this is the case for the Gaufi constraint operator. However, often the physically interesting 
operators are unbounded and in this case they are only deflned on a dense subspace of 
"Hkin- Furthermore, in general, the spectrum of the operators will have a continuous part 
and the associated (generalized) eigenvectors will no longer belong to Hkin- In our case 
this situation occurs whenever zero will not lie in the discrete part of the spectrum of the 
constraint operators. A similar situation occurs already in standard quantum mechanics 
for the momentum and position operators both being unbounded. The spectrum of the 
momentum operator p = —ih-j^ is the real line and the (generalized) eigenfunctions are 
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plane waves -^/t {x) = e**^^ that are not square integrable on (M, dx) and thus no element 
of T-Lqm = i2(R,rf2:)- Likewise, the generalized eigenfunctions of the position operator 
are delta functions and therefore also no elements oH-Lqm- In those cases one has to look 
for solutions in a larger space than 'Hqm and, in the case of LQG, ?^kin respectively. The 
mathematical framework that can be used here are so called rigged Hilbert spaces (also 
called Gelfand triples). They consists of a sequence of spaces 

VcncV (101) 

where I? is a dense subspace of Ji endowed with its own intrinsic topology, that is assumed 
to be stronger than the one induced from T-l and 2?^ is the dual space of D containing 
all continuous anti-linear functionals on V. The dense subspace is usually chosen by 
the requirement that it is the largest dense subspace of H such that it is an invariant 
domain for arbitrary powers of the elementary operators. For instance in the case of 
quantum mechanics the dense subspace is the 5(K), the Schwartz space of smooth rapidly 
decreasing functions on R. Associated with a given Gelfand triple there always exist a 
second rigged Hilbert space 

VancV' (102) 

which consists of P,?^ and the dual space V of all linear and continuous functionals 
on V and it is the second rigged Hilbert space that we wiU discuss for LQG. In the 
case of quantum mechanics when the spectrum has a continuous part the associated 
(generalized) eigenvectors expressed as Dirac ket vectors live in I?^, whereas the home of 
the corresponding Dirac bra vectors is the space V . For LQG we restrict our discussion 
to the latter Gelfand' triplet in (J102I) . Let us denote with 2?kin the dense subspace of "Hkin 
and as we saw before elements of I?kin are cylindrical functions. Apart from the dense 
subspace I?kin itself the rigged Hilbert space framework requires that 2?kin is endowed 
with an intrinsic topology stronger than the one induced from 'Hkin and as far as physics 
is concerned it does not yield a particular choice of a topology here. In order to avoid such 
a choice at this stage in LQG we do not consider the topological dual space T)'^^^ but the 
algebraic dual instead denoted by 2?ki,j. Then we look for solutions in the algebraic dual 
■^kin that is the space of all linear but not necessarily continuous functionals £ on I?kuj 



■^kin ^^ naturally equipped with the weak *-topology of pointwise convergence of net^ 
As before we have the following topological inclusion 

2?kin C Hkin C V^,^ (103) 

since any functional that converges strongly in the norm of 'Hkin will also converge point- 
wise. In order to formulate a requirement for solutions to the constraints in 2?kj„ we need 
to extend the action of the operators from Hkin onto T)^^^. On those linear functionals £ 
that lie in Hkin we want the action of the dual operator to agree with the usual one on 
"Hkin. Let us denote the extension of an operator O by O' then we define 

&£] (/) := £(6V) (104) 

where f denotes the adjoint in 'Hkin- Now suppose we have an £ e 'Hkin C Pkin ^^'^n using 
Riesz representation theorem we find a unique fg G Hkin such that i can be expressed as 
£ = {fg, .)>^^i„ where (., .)Hki„ denotes the inner product on ?^kin- Then we obtain 

d'£] {f)^£{dU) = {h,6^fU,^„ = {d.fiJ)n,,^, (105) 

which explains the use of the adjoint in equation (|104p . Looking for solutions in 1)^^^^ 
corresponds to finding linear functionals that satisfy the following requirement 



C"£ 



(/)=£((7t/)^0 for all C, / e 2?kin (106) 



A net, that is a generalization of a sequence, {£") converges in V^^^ if {£a{f)) converges to £{f) for any / 
in ©kin- 
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Let us denote the space of solutions to the constraints by 2?*)^ . In general physical oper- 
ators, those that commute with all constraint operators, will be unbounded and therefore 
only be defined on a dense subspace of the physical Hilbert space "Hphys denoted by Pphys, 
which is an invariant domain for the algebra of physical operators. The set T)*^ is 
then the algebraic dual of I'phys and as before again one has also at the physical level a 
topological inclusion 

^phys C "Hphys C 2?phys (107) 

Given the situation that we have found I^p^yg there exists a systematic way to construct 
an inner product by means of a so called rigging map 77 

r, : Pkin ^ 2?;hys, / ^ Vif) (108) 

that maps elements in Pkin hito 'D*^^ . ri{f) being an element of 2?*j is a linear functional 
on 2?kin that additionally satisfies the condition in equation (|105p . hence [rj{f)]{f) is a 
complex number and if in addition 



(*) WW) = wm) (109) 

i^^) [^(/)](/)>0 for ah/, Wm) = ^ f = (110) 

(til) [vifMc^g + ^h) = a^mg) + P[rj{f)]{h) (111) 

are satisfied we can use [??(/)](/) do define a physical inner product on the image of rj 
given by ?7(I?kin) G ^phys- The bar in (i) denotes complex conjugation. Condition (iii) is 
trivially satisfied since [?](/)] is a linear functional. Whether (i, ii) are satisfied and thus 
[77(.)](.) defines a positive semi definite sesquilinear form needs to be checked once the 
explicit form of rj has been constructed. Apart from the minimal requirements (i-iii) we 
want the inner product to satisfy the property that adjoints with respect to the physical 
inner product (., .)phys represent the adjoints in the corresponding kinematical case. That 
means that we can either first extend the operators to T)*^^ and then construct the adjoint 
or take the adjoint first in Pkin and then extend the adjoint operator to T)*^ and we will 
obtain the same result. This yields the following condition for the physical inner product 

(IV) ((6')V,Vi)phys=((0^)V>)phys (112) 

Here we denoted the adjoint with respect to the physical inner product with the same 
symbol f we used in the kinematical case. Furthermore, the rigging map needs to be 
constructed in such a way that the physical operators O defined on Pphys respectively 
preserve the space of solutions X'pWs, that means 

{v) d'v{f)^7j{df) for all / e Pkin (113) 

A physical inner product on the image of ?/, that is ?/(I?kin) C T)^^^ can then be defined as 

{^, v;)phys - {vif),vif))pbys ■■= mm) (114) 

The physical Hilbert space Hphys is then defined as the completion of rj with respect to 
(•J •)phys With that definition of an inner product the condition (iv) in (|112p is automati- 
cally satisfied 

((6')^V'>)phys = (^,OV)phys 

= {v{.f),V{Of))pl.ys 

= iviomf) = [o'viMf) = iviMo^f) 

= {v{0^f),Vif)Uys 

= ((Ot)V,^)phy. (115) 
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An example how a rigging map can be explicitly constructed is the so called group aver- 
aging procedure. Consider a set of constraint operators C/ where / labels the individual 
constraints that are self-adjoint and their first class algebra is a Lie algebra, then we can 
define unitary operators by using the exponential map 

U{g)=expLY,(^'Ci] (116) 

with parameters 0^ g T C M yielding a unitary representation of the Lie group. The 
condition that a function / is annihilated by the infinitesimal generators of the constraints 
carries over to the requirement that a linear functional is a solution of the constraints if 
it is invariant under the action of those unitary operators, that is for / £ I?kin 

mi9)Y]if)^m{9)yf)^m for aU .gee (117) 

Using the unitary operators we can define a projector on physical states given by 

P= fd^iH{g)U{g) (118) 



where ^h denotes Haar measure on G. The property of the Haar measure that it is 
invariant under right and left translations are important for showing that the projector 
defined above indeed projects on physical states as can be seen below 

Ui9)Pf - U{g) J dtiH{h)U{h)f 

G 

dfiH{h)U{g)U{h)f 

G 

d^iH{h)U{gh)f 

G 



dfiH{g-'~g)Ui~g)f 

G 

= Pf (119) 

In the third line we used that G is a group, in the fourth line we introduced the new 
integration variable g := gh and in the one before the last line we used that fj,h is invariant 
under translations. The rigging map rj can now be expressed in terms of the projector 
defined in pTS)) 

ry : I?ki„ ^ Hphys C 2?;hys f^v{f)-^Jdfi{g)H{U{g)f,.)Mn (120) 

G 

and the physical inner product can then be defined as 

('7(/),^(/))phys:=b7(/)](/) (121) 

Likewise to the case of quantum mechanics where the Dirac bra and ket vectors in the 
spaces V^ and V are distributions on V also here the linear functional ?y(/) defined 
above is a distribution on I?kin- In this sense the rigged Hilbert space framework com- 
bines Hilbert spaces with distribution theory and allows us to understand the case of 
unbounded operators with continuous spectra along the lines of the standard language 
used in quantum theory. 
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Before we will apply RAQ to LQG in order to solve the diffeomorphism constraint, we 
will discuss a simple example from quantum mechanics, where one basically could guess 
the physical Hilbert space immediately so that the individual steps of the RAQ program 
can be understood in a simple situation. We consider a two dimensional quantum me- 
chanical system with kinematical Hilbert space 'Hkin = L2{M.'^ , dx"^) with one constraint 
given hy C = pi. In the first step of RAQ we need to choose a dense subspace P^jn of 
^kin on which the elementary operators Xi,X2,Pi,P2 and arbitrary powers of them can 
be defined. As usual we choose 5(K^), the Schwartz space of rapidly decreasing functions 
on R^ for 2?kin- Now we are looking for solutions to the constraint C = pi = —ihd/dxi. 
These are functions that do not depend on the variable xi and since the inner product of 
"Hkin involves an integral over R of dxi those functions are not normalizable and thus no 
elements of Hkin- However, they can be mathematically defined as elements of 2?ki„ and 
hence linear functionals £ on 2?kin defined as 



£ : Pkin ^ C / ^ £{f) := {£, />km - / d'x£ix)f{x) (122) 

The second step of RAQ consists of finding solutions in V^^^ and in our case these are 
functionals that satisfy 

[C'£]{f) := £{C^f) ^ J d:'x£{x)p{f{x) = th J dxidx2-^^{xi,X2)f{xi,X2) = (123) 

R2 R2 

Obviously, the elements of 'D*^^ C Pkin are those linear functionals that do not depend on 
xi. Therefore, physical operators defined on Hphys or a dense subspace X'phys respectively 
are operators that depend only on X2,Pi,P2 because those commute with the constraint 
C = pi. Now, since pi is a constraint operator involving pi correspond to classical objects 
that vanish on the constraint surface and are thus rather uninteresting physical operators. 
For this reason the relevant physical operators will be functions of X2,P2 only. Hence, the 
physical Hilbert space in this simple example is just "Hphys = ^2(1^1 dx2)- Let us show that 
we end up with the same result when we use the rigging map to construct the physical 
Hilbert space. The projector for group averaging in this example has the form 

P :./..* (124) 

and using P to define the corresponding rigging map r/ : I?kin — ^ ^phys ^^ obtain 



7y(/) := / dt{e'*P'f,.)k^n ^ dt d'xf{xi+t,X2) (125) 



where we used that the unitary operator exp{itpi) generates translations. The physical 
inner product can then be defined as 



('7(/), ^(/))phys := Wm)^ J dtj dxidX2f{xi+t,X2).f{xi,X2) (126) 

K R2 

Let us introduce the variable t' := t + xi then we rewrite the expression in (jl26p as 

ivif), Vif))pl.ys ^ JdX2 ijdt'fit',X2) ljdx,fixi,X2)] = (Vi, ^)L.(R,d..) (127) 
R \r / \R / 

showing that the physical Hilbert space we end up with is indeed Hphys == ^2(1^! dx2). 
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Finally, let us say a few words about alternative methods to solve the constraints in 
cases where the refined algebraic quantization (RAQ) framework cannot be applied. This 
is for instance the case for the Hamiltonian constraints since their algebra involves struc- 
ture functions instead of structure constants. A way to circumvent this problem is to 
replace the Hamiltonian constraint by the corresponding master constraint M defined in 
(ED) |27, , that (i) satisfies a trivial algebra and is constructed in such a way that one 
obtains a self-adjoint operator [ST]. Now one could on the one hand apply RAQ to find 
solutions to the Master constraint. On the other hand one can also use the direct integral 
decomposition (DID) method in order to construct the physical Hilbert space. In the 
latter case one uses the fact that for self adjoint operators O on separable Hilbert spaces 
the Hilbert space can be written as a direct integral 



U ~ / dfi{\)Hx (128) 



whereby fj, is the spectral measure and "Ha are again separable Hilbert spaces with an 
inner product induced from H. On each of this Hilbert spaces the operator O acts by 
multiplication. So assuming that O is a constraint operator the physical Hilbert space 
is just the Hilbert space Hx associated with the eigenvalue A = 0, that is Hphys — ^o- 
For LQG the kinematical Hilbert space on which the master constraint operator M will 
be defined is non separable. However, T^kin decomposes into a direct sum of separable 
Hilbert spaces each of which are left invariant by the action of M and therefore the DID 
method can be applied to each of the individual Hilbert spaces in the direct sum, see [52] 
for an application of this method. 

3.3.3 Quantization and solutions of the diffeomorphism constraint 

When looking for solutions of the diffeomorphism constraint, we would like to proceed 
in a similar way as for the Gauss constraint. Doing so, the first step consists in defining 
the infinitesimal version of the spatial diffeomorphism operator which is classically of the 
form 

C{N)^ jd\N^Fi,E]{x) (129) 

s 

where N'^ is the shift vector. Fab — Flb'^j/'^ the curvature associated to A^ and E"^ the 
densitized triad, where we neglected the terms involving the Gauss constraint Gj. Consid- 
ering the term involving the Gauss constraint it is possible to rewrite C{N) as an integral 
over S oiE°i and {Cj^A^)a yielding, similar to the Gauss constraint, to an expression that 
involves the densitized triad smeared over a 3 dimensional integral with a phase space 
dependent smearing function. Now, when one tries to proceed and expresses C{N) in 
terms of holonomies and fluxes one realizes that the operator corresponding to C{N) does 
not exist on T^kin- Note that recently in [53] using the habitat introduced in [75] a quan- 
tization of the infinitesimal diffeomorphism operators was performed using a densitized 
shift vector of the form N'^ = E^ /^ yielding an expression that is less singular and can 
be quantized. As we will see in section |4?2] the same is also true for the classical expression 

q'^'CaCb. 

However, that Ca itself does not exist is not a problem as far as the Dirac quantization 
program is concerned because the requirement that the representation should be diffeo- 
morphism covariant ensures that finite spatial diffeomorphisms are implemented unitarily 
(see also section \3A\ for a more detailed discussion). As a consequence we can work with 
the operators denoted by U{f) for ip G Diff(I]) instead when looking for solutions of the 
diffeomorphism constraint. 
In order to define the action of U{ip) on a given spin network function we define the 
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multilabel s — {7, j, I}. Then we have for Ts = T -> f 

Ui^)Ts = T^.s f.a. ip e i?^//(S) (130) 

with ip-s — {ip-j,ip-j,ip-I} and 

V ^'- J eeEi-yis)) 

ifi-I := ( ((y9 • /)^(„) = ^ ) (131) 

\ / veV('y{s)) 

where E{'-f{s)) denotes the set of edges of the underlymg graph 7(3) of Tg and ^(7(5)) 
denotes the set of vertices of 7(s). We see that the action of U{(p) maps the graph 7(3) 
on which Tg is defined onto its image under the diffeomorphism (p. Note that definitions 
above are understood in the way that the order in which the edges at the vertices are 
coupled to obtain the corresponding intertwiners is not changed under the action of the 
diffeomorphism operator. That the infinitesimal constraints do not exist as operators can 
be also seen from the fact that the family of operators U{ip) are not weakly continuous. 
A family of unitary operators U{ipt) is said to be weakly continuous if we have that 

\im{Ts,U{cpt)Ts,)^{Ts,Tg,) f.a. T,,T,, (.U^^ (132) 

Let tfY be a one-parameter family of diffeomorphisms generated by a vector field V ^ Q. 
If we choose 7 in the support of V then there exists e > such that for all t G (0, e) 

ft il) 7^ 7- Now we choose Ts ~ Tg' and obtain 

lim(T, , U{^Y)Ts) - lim(T, , T^v^,)) = ^ 1 = (T, , T,) (133) 

which shows that the finite diffeomorphisms are not weakly continuous. For the reason 
that any cylindrical function /^ can be written as a linear combination of spin network 
functions, the action of U{^) can be easily extended to /^. In contrast to the Gauss 
constraint, solutions to the spatial diffeomorphism constraint will not live in a subspace 
of 'Hkin and we need to apply the RAQ program here, which has been done in [H]. An 
application to 2+1 Euclidian gravity can be found in [Bl] and in [5S] the RAQ quantization 
program was applied to a scalar field toy model. 

As a first step in RAQ we need to choose a dense subspace of Hkin on which arbitrary 
finite powers of the elementary operators are defined. In the case of LQG I?kin = Cyl the 
space of smooth cylindrical functions. Now we are looking for solutions in the algebraic 
dual 2?^ijj = Cyl* and for this purpose the action of operators and particularly of U{(p) 
on ?^kin needs to be extended to functionals £ in 2?kin ^^'^ i^ given by 

[U'iip)e]{f)^iiU-Hf)f) f.a. /ePki„,¥'e^*//(S) (134) 

where we denoted the extension of the operator with a prime as above. This allows us to 
formulate the requirement for solutions £ 

[U'{^)£]{f)^l{U-\^)f)l£{f) f.a. feV^,^,^eDiff{T.) (135) 

Those solutions are elements of 2?phys C PJijj defined as the set of those £ G P^jj^ that 
satisfy the condition in (|135p . Since the spin network functions lie dense in I?kin = Cyl we 
can restrict the construction of solutions to the diffeomorphism constraint to them and 
then express diffeomorphism invariant cylindrical functions as linear combinations of the 
solutions associated with spin network functions. 

[t/'(^)^](r ^j) - t{U-\ip)T^ jf) = ^(T^-i(^) J /) = i{Ts) f.a. / e Pkin, f G DiffiE) 

(136) 
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Let us introduce the orbits of s associated with Dif f{Yi) denoted by [s] 

[s]:={vp-s, 0ei?z//(E)} (137) 

For aU orbits [s\ a diffcomorphism invariant distribution can then be constructed as 

ks]= E^^''''-)^- (138) 

s'e[s\ 

We need to ensure that the sum in (|138p converges when applied to some spin network 
function T, as otherwise ll[s] would not be an element of T)*-^ . Fortunately, due to the 
orthogonality property of spin network functions that follows from the Peter & Weyl 
theorem, already mentioned in section [3. 1[ in the expression 

ks] (Ts) = Y. <^-' ' ^^>ki„ = X[s] i~s) (139) 

s'e[s] 

only one term will be non- vanishing where x denotes the characteristic function. First, 
we only obtain a nonvanishing contribution if s € [s] , meaning the graph T^ is defined on 
must be diffeomorphic to the graph associated with the multi-label s. Furthermore the 
righthand side of (|139p is only nonvanishing if additionally the representations associated 
with the edges agree. Note that in the context of group averaging (|139p can be understood 
as group averaging over the orbits of [s] under Diff(E) using a counting measure on the 
orbit. That those £[5] are invariant under the action U{ip) can be easily seen 

u'{ip)i[4Ts) = E(r,.,f7t(^)r,)ki„ 

s'e[s] 

= E(f7Mr,.,r,)ki„ 

s'e[s] 

= 2^ {Ts", Ts)kin 
s"e[s] 

= ks]iTs) (140) 

As far as the operators corresponding to strong Dirac observables with respect to spatial 
diffeomorphism are concerned, these are operators that strongly commute with U{ip), one 
can show 44J that "Hkin splits into mutually orthogonal super selection sectors, that is, it 
decomposes into a direct sum of Hilbert spaces associated with the individual orbits [s] 

Hki„ = 0^|;l with n\:i=^nL (wi) 

and (strongly) diffeomorphism invariant operators preserve those individual Hilbert spaces 

[si [si 

^kin- The rigging map rj is then constructed for each individual "Hj^^ separately. For this 

[si [si 

purpose let us choose some dense subspace I'kjn C Hki„ and consider the algebraic dual 

r 1 * r 1 r 1 * 

(■^kin) ■ ^ family of rigging maps i]a ■ 2?ki,j — >■ (2?jjff ) can then be defined as 

?ya(r,) :=«[,]£[,] (T,) (142) 

where au] > G M is some up to now unspecified number labeling the family of maps. 
However, in general aui depends on the orbit of the multilabel s and not on the orbit of 
the graph 7(3) only. As a consequence we will not be able to define a symmetric inner 
product by means of rja unless we will choose the coefficients a^g] in such a way that they 
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depend on the equivalence class of the graphs only. 

In order to discuss how the coefficients ar^i need to be modified and also make the connec- 
tion to the group averaging more transparent we will rewrite rja as sum over diffeomor- 
phisms in DifF(E). Note that the naive ansatz where one averages over the total group 
Difi^(I]) does not work for the reason that there exist uncountably infinitely many diffeo- 
morphisms that leave a given spin network T^ invariant. Thus the trivial difFeomorphisms 

need to be factored out in order to obtain a well defined element of (2?jjff )* C (I'jjin) • We 
consider the diffeomorphisms that leave a given spin network function Ts invariant and 
distinguish between two cases; 

TDi//[,] := {ipeDiffiJ:)\ip-s^s} 

TDiff^^] := {ipeDiffiJ:)\^-j^j, ^ie) = e f.a.eeEi^is))} (143) 

where ip{e) — e should be interpreted as saying that the edge including its orientation 
are preserved. We used a notation analogous to the one introduced in [S]. However, in 
[5] the group averaging and the rigging map is defined in terms of equivalence classes of 
graphs and not of spin network functions. Both TDiff^g] and TDiff^^] are subgroups 
of Diff{T,). TDiff^s] includes all diffeomorphisms that preserve the spin network Tg 
whereas elements of TDif fhy-i additionally preserve all edges of the underlying graph 7(5) 
and hence do not depend on the multilabel s anymore but only on the graph. The quotient 
of the two groups turns out to be a finite group, the group of spin network symmetries of 
Ts denoted by 

S[s] ■■= TDiff[s]/TDiffi^] (144) 

The group averaging can then be expressed as an averaging over the group of spin network 
symmetries and an averaging over the remaining difFeomorphisms that fill the orbit [s] of 
s which is Diff{'S)/TDiff[s] — Z?i//[s](S). The projector associated to the averaging 
over Su] can be expressed as 



\Su 



i^eSi, 



Note, that in general the size of S[s\ and S^s'] will be be different even when the graphs 
associated with the multilabels s and s' are diffeomorphic. As a consequence, we need 
to consider this in the definition of the prefactor that enters into the rigging map rja to 
ensure that the prefactor depends on the equivalence class of the graphs only. Hence, 
when we consider also the averaging with respect to those difFeomorphisms that move the 
spin network we obtain for the rigging map 

i^aiTs) ^ a^-y]\S^s]\ E (f>M^S[„r,,.)ki„ (146) 

<P6-D*//[,1(S) 

Here we introduced a^^j as the real number that one obtains from defining the group av- 
eraging in terms of graphs and not spin network functions in order to make a comparison 
with the notation used in [5j possible. We have defined the rigging map above in a way 
that the norm of the most unsymmetric spin network function is given by a\^^ which can 
be chosen to be one. 

A diffeomorphism invariant inner product can the be constructed on the image of rja 
denoted by '7a(2?{;i„) C {T^ll^) and the diffinvariant Hilbert space HJ^jg is then the com- 
pletion of 2?|iifj with respect to (., .)diff 

(r;,(T,0,rya(T,))diff := [Tja{Ts)]{Ts-) = a[^]|5[,]| ^ (C/(^)Ps,.jT, , T.^kin (147) 
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It can be shown that the so defined inner product satisfies ah five requirements mentioned 
in section [5.3.21 in equations (|109p . (|112p and (|113|) . The Hilbert space Hdiff can then be 
constructed as the direct sum of the individual Hilbert spaces T-L^iS- Those requirements 
in (|109p . (|112[) and (|113[) do not yield any further restrictions on the factor a[^] in (|147p ex- 
cept the already implemented condition that a[^] is a real positive number. Consequently, 
as far as operators are considered that correspond to strong Dirac observables with respect 
to spatial diffeomorphisms the inner product of the solution space is not unique because 
we have a freedom to normalize the inner products associated with the mutually orthogo- 
nal spaces T-L^^ff- However, one expects that operators corresponding to weak observables, 
that is, operators that commute with the spatial diffeomorphism constraint only modulo 
constraint operators, will mix super selection sectors and could therefore yield additional 
requirements for the normalization constants and hence further restrict the ambiguity of 
the inner product on 'Hdiff. 

Finally, let us briefiy mention that also the Gauss constraint can be solved using group 
averaging. Let us denote the unitary operators corresponding to finite SU(2) gauge trans- 
formations by U{g), then the associated rigging map is given by 

v{Ts)= J l[d^^H{{9(x)){U{g)T,,.)k^n (148) 

Using T] above to construct the gauge invariant Hilbert space "H^j^ yields to exactly the 
same result as we obtained in section 13.3.11 Note, that solving the Gauss constraint can 
equivalently be done either before solving the diffeomorphism constraints or afterwards. 
The resulting Hilbert space "H^ijj is, similar to the kinematical Hilbert space 'Hkin still 
non-separable. In _B6^ a proposal was introduced to obtain a separable Hilbert space for 
'Hdiff by allowing fields to have isolated point of non-differentiability. The associated gauge 
symmetry is an extension of Diff(I]) and includes homeomorphisms that, together with 
their inverse are smooth everywhere except at at most finitely many points. An application 
of the RAQ program with respect to this extension of Diff(I]) yields a separable HdiS- 

3.4 More on representations of 21 

We have already discussed one representation of the canonical commutation relations 
(|4T|) (or rather their integrated version ((53|) ). This was the diffeomorphisms covariant 
representation (JSQ]) on the Ashtekar-Lewandowski Hilbert space Hkin- In this section we 
would like to discuss other representations. Why would it be interesting to do so? There 
are some aspects of the standard representation that are peculiar, and one might wonder 
whether there are other representations that do not have these properties. Examples are: 
In the standard representation, 

• there are no operators representing the connection A directly, 

• the generator of the unitary maps implementing diffeomorphisms does not exist, 

• the spectrum of geometric operators is purely discrete. 
Natural requirements for a "fundamental" representation include: 

1. The representation is irreducible, i.e., each state in the representation Hilbert space 
is cychc. 

2. The representation is diffeomorphism covariant, i.e., there exist unitary operators 
implementing the diffeomorphisms. 

3. There is at least one state that is invariant under diffeomorphisms. 

A representation with all these characteristics is, however, equivalent to the the Ashtekar- 
Lewandowski representation (j59p . This is a consequence of the uniqueness theorem 
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Proposition 1 (|67[ 145] ). There is only one cyclic representation of ^ with diffeomor- 
phism invariant cyclic vector - the Ashtekar-Lewandowski representation ([5^ 



We note that to really prove the above proposition, the algebra 21 and the class of 
diffeomorphisms has to be defined in great detail. We also note that there are interesting 
representations that violate one or more premises of the above proposition|3 For example, 
[68] exposes a difFeomorphism invariant representation that is reducible. In j69], a cyclic 
and diffeomorphism invariant representatiopn of an algebra slightly larger than the one 
from the above proposition are described. 

In the following, we will however describe a simpler class of representations due to 
Koslowski, [70l 1711172] . These representations have been called representations with classi- 
cal spatial background geometry. The basic observation is that the Ashtekar-Lewandowski 
representation is a member of a family of very similar representations. The AL ground 
state is peaked on a degenerate triad E^^") = 0. But it is easy to construct similar rep- 
resentations which, however are peaked on different classical triad fields E^^' . A rough 
sketch of the situation is as follows: 

AL ground state ft = So{E) in momentum rep. (149) 

Ground states w. background flEW = S^(o) (E) in momentum rep. (150) 

How are these representations defined precisely? The Hilbert space, and the representation 
of the holonomies stays the same, 

■Hb(o) = "Hkin, T^Ewih) = h, (151) 

but the representation of the fluxes is changed by adding a c-number term, 

^£(0, {Er,{S)) = Xr,{S) + E^^\S) id, El^\S) = I *Ef^n' d\. (152) 



JS 

The new representations have the following properties. 

• The representations 7r^(o) are cyclic. 

• One can regularize the geometric operators in the representations 7r^(o) with exactly 
the same methods as in the standard representation (described in sec. 13. 2p . The 
resulting operators are 

MS)^Arp,i^{S) + Ar{E^°\S)\d, V [R) = Vp,i^{R) + V {E^°\ R) id, (153) 

where 5 is a surface, R a region, the subscript AL denotes operator in the standard 
representatiorlj, and the c-number terms are given by the classical value in the 
respective background. This result is very simple to state, but the proof is non- 
trivial [TT] . 

• Only the symmetries of E^-^' can be implemented unitarily in the new representa- 
tions. 

We also note that for i?*-°^ = 0, we recover the standard representation from ()151|1 . ()152|) . 
Equation (|153p shows that the new representations can be interpreted as containing a 
"geometric condensate" . 

If we want to use the more general representations in the place of the standard rep- 
resentation when implementing constraints, the failure of the diffeomorphisms to be im- 
plemented unitarily is of concern. For example, for the operators C/((/') implementing the 
diffeomorphisms in the standard representation one finds that generically 

C/((/.)7rs(o) (S„(5))t/(0)t ^ ^s<o, {E^,n{<i>{S))) (154) 



®To be precise, the examples given in the literature are for the somewhat simpler case of the structure group 
SU(2) replaced by U(l), but it is likely that they generalize to the SU(2) case. 

Val can be either, the Ashtekar-Lewandowski, or the Rovelli-Smolin version, depending on which regular- 
ization procedure is chosen. 
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for difFeomorphisms (f>. But one can easily show that one can also not find other unitaries 
that will do the job. The reason is that the geometrical background iS^^^ is fixed can 
can not transform under any operation on T-Le^") — ^kin- There is a way to remedy 
this problem, at the price of going over to a much larger Hilbert space and a reducible 
representation. Let us use the notation 

|T)s<o, ^|r,i?(°)) (155) 

for a spin network function T interpreted as an element of Hew ■ Then define 

E(0)e[-E<0)] 

Here the direct sum is over all background triads in one gauge and diffeomorphism equiv- 
alence class, 

:E(oTe[^(o)] ^=^ 3g,0: ;bW =Adg((^*^(°)), (157) 

where g denotes a gauge transformation and (j) a diffeomorphism. We note that the Hilbert 
space ()156p is thus labeled by a spatial metric modulo diffeomormphisms, that is, a spatial 
geometry, or a point in superspace. 2t is represented on T-LiEiO)] through the direct sum of 

the representations tt-^^^ for i?(o) e [£^^"^1, and gauge transformations, diffeomorphisms, 
as well as their semidirect product can be implemented unitarily. The representation on 
H^EW] is not cyclic, however. 

Starting from 'Hie(o)]^ it is possible to obtain states that solve the diffeomorphism and 
Gauss constraint by group averaging, just as in the case of the standard representation 
discussed in detail in sections 13.3.21 Let us consider the diffeomorphisms as an example. 
We make the definitions 

TDiff([3],£;(o)) := {Diffeos (/) : (j)^E^°^ = E^^^ and • s = s} (158) 

TDiff([^]^^(o)) := {Diffeos : </.*^(") = ^(°) and 0(e) = e for ah edges e of 7} (159) 

Note that these definitions exactly parallel those used in the detailed treatment of the 
group averaging procedure in section 13.3.31 Like in the standard case both of these sets 
actually form groups. Moreover their quotient , the group of spin network symmetries 
given by 

S([s],_B(o)) = TDiff([^]^^(o))/TDiff(^^^(o)) (160) 

can be shown to be a finite group. Here, it is understood that Ts is cylindrical on the 
graph 7. Then group averaging for a spin net T^ in the iJ^^^-sector of HtEW^ effectively 
reduces to averaging over diffeomorphisms modulo TDiff/r^i^^(o)-) denoted by Diff[s](E) = 
Diff(E)/TDiff/r^i^^(o)') and over S(-ui £;(o)). More precisely, let 

(T„i?(0)|T,-,i?'(0)) :=a[,]|5M| ^ {T,,E^^) \Ps^^^^^^,,Ul\T,,,E'^''^), (161) 

G(0)GDiff[,](S) 

where the projection Pq ,„> is defined as 

(1s1,e('')) 

Ps , |r„£;W):= ^ Y. U^Ts^E^'^^). (162) 

Now it is easy to show that 

Proposition 2. The linear functionals (T,,-^^^'! are well defined, finite, and diffeomor- 
phism invariant, 

(T„ ^(°) I o C/^ = (r„ S(o) I for all <i> G Diff. (163) 
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Again, these definitions reduce, for the case i?*-''' = 0, to the diffeomorphism averaging 
map defined in section 13.3.31 Similar resuhs can be obtained for gauge transformations 
and, taking the semidirect product of diffeomorphisms and gauge transformations, for 
bundle automorphisms. This means that the quantum kinematics can be developed to 
the same point for the new representations as for the standard representation. Among 
other things, this shows nicely that diffeomorphism invariance is not antithetical to being 
peaked on a fixed geometry. 

As an example consider the operator Vol for the volume of the entire spatial slice. 
Also in the new representations, it commutes with all automorphisms. It thus defines an 
operator on 'Haut- Moreover, this operator acts in precisely the way one would expect. If 
/ is an eigenstate of the volume operator in the standard representation, with eigenvalue 
A, then \f, E''^"^) is an eigenvector of Vol with eigenvalue A + Vol'^^\ 

We also note that the invariant states resulting from the vacuum representation are 
again contained in these constructions as the special case £''"' = 0. Finally, while -ffrsto)] 
is large, this is partially remedied by the group averaging. For example, vectors | 1,E^^'), 
1 1,E''-^') S Hi^(of\ are mapped onto the same vector in "Haut- More generally 

(/,^.i?(°)| = (CV,i^^°^|. (164) 

4 Quantum dynamics 

In this section we discuss the quantum dynamics of loop quantum gravity. In the case that 
we use pure Dirac quantization the dynamics is encoded in the Hamiltonian constraint 
and its quantization is explained in section 14.11 When we consider the Brown-Kuchar or 
the scalar field model earlier introduced the dynamics is encoded in a so called physical 
Hamiltonian. However, the quantization of the latter relies in both models on techniques 
that have been used for the Hamiltonian constraint. The details of the quantization of 
the Brown-Kuchar model are discussed in section [4. 21 whereas the quantum theory of the 
scalar field model is presented in 14.31 



4.1 The quantum Hamilton constraint 

As we have seen in section [2?T1 the Hamilton constraint of the classical theory is given by 



^^EfE^e^hF^,-P{l + l3')^. 



C = ^ —E^E^e^^kF^, -/3(1 + /3') —E^E^K^ . (165) 



= :Ce =:T 

In the present section we will discuss how to turn this classical expression into a well defined 
operator. The general difficulty with this is obviously that C is a complicated nonlinear 
function in the phase space variables, hence ordering problems present themselves. There 
are also some specific difficulties with the expression: 

• (I165P contains the inverse volume element. The volume element itself has a large 
kernel when quantized, see the discussion in section [3.21 so its inverse is ill defined. 

• The expression (jl65p contains the curvature F oi A, as well as the extrinsic curvature 
K. For neither of them there is a simple operator in the quantum theory. 

A guiding principle in the quantization process can be the Dirac algebra p8)) - ([20)) . In par- 
ticular, the quantum Hamiltonian constraint should be invariant under gauge transforma- 
tions, covariant under diffeomorphisms, and the commutator of two Hamilton constraints 
should give a diffeomorphism constraint. 

We should say that the knowledge about the quantization and implementation of the 
Hamilton constraint is not complete. Many things remain to be understood. But we will 
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show that at least there is a strategy that leads to well defined constraint operators. Given 
the difficulties outlined above this is highly nontrivial in itself. 

The quantization strategy we will describe in the following is due to Thiemann |73l 
[74l [75] , but draws on important earlier work and ideas by Rovelli, Smolin, Lewandowski 
and others. Our presentation is in part based on [S]. 

Thiemann's tricks 

The quantization is based on two key ideas. The first one is to use various ingenious 
classical identities to express parts of the Hamilton constraint in terms of Poisson brackets 
before quantization. The second one is to express curvature in terms of holonomies. Let 
us explain them in turn. 
Let 

V 

It, Jt, 

be the total volume of the spatial slice, and the integrated extrinsic curvature. Then 



( d^x v/det q, K=fd^xKlE^ (166) 



' ^ ^_e-''^{t/,A,n, Ki = -{K,Ai}. (167) 



Vdet q K K 

These identities can be used to write 



CEiN) = c / d-'x iVe-^"^ tr {Fab{A„ V}) , (168) 

T{N)^c' f d''xNe'''>-tY{{Aa,K}{At,K}{A,,K}), (169) 

where we have used the notation for the two parts of the Hamilton constraint introduced in 
(jl65|) . The constants are c = 8/k and c' = 16/k^. The idea behind these reformulations is 
that it is natural to replace Poisson brackets by commutators in the quantization process, 

{•,-}~-^[-,-]. (170) 

This means that the quantization would be greatly simplified if operators existed for the 
quantities V, K. Indeed we have already seen in section [X2] that an operator exists for V. 
With respect to K, the identity 

K = {V,Ce} (171) 

suggests to first quantize Ce, and then use the commutator with the volume operator 
to define the operator for K. Thus we have already dealt with two of the difficulties 
regarding the quantization of C: The inverse volume element is gone, and the extrinsic 
curvature is dealt with. What remains is the quantization of the curvature F of A. Here 
we use the well known fact that holonomies encode information about curvature. Let S 
be an oriented surface such that the integral L F is small, and let a be the (oriented) 
boundary of S. Then the first term on the right hand side of 

^^ = i(v-/.;Vo((//)') (172) 

is a good approximation to the left hand side. Let e be an edge starting at a point s(e). 
A similar approximation plus a second Taylor expansion gives 

ee-{s{e)){A{s{e))a,V} « h-'{K, V} (173) 

where e is the tangent to e in a chosen parametrization e{t) , and e is the coordinate length 
e = J dt oi the edge in the given parametrization. In this way, we can express curvatures 
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and connections by holonomies. Putting everything together we can get a Riemann sum 
approximation of the Euchdean part of the constraint, 

3 

CE(iV) « 4°) := ^E^(«n)E [(C(n) - V(n)) /^;/(n) {hs,in),V}] . (174) 
n 1=1 

Here {D} is a decomposition of S into 3-dimensional cells, and for each ceh a point v\j has 
been fixed. {Q;/(n)} is a set of loops and {s/(n)} a set of edges such that their tangents 
span the tangent space in the point va in the following sense: there is a basis {^/(n)} 
of the tangent space at !>(□), such that bj{C\) is tangent to both ai{0) and S7(n), and 
compatible with their orientations. We call the data ({□}, {vn}, {•?/(□)}, {^/(n)}) a 
regulator of Ce, and sometimes denote it simply by D. The exact shape of these cells, 
loops and edges does not matter. The approximation is good as long as the cells are much 
smaller than the scale on which the fields A, E vary, and the loops and edges stay within 
the cell. 

Finally, we can consider families of regulators such that the cells shrink to points. Then 
the corresponding approximations will converge to the exact result for a wide variety of 
such families. 

The same kind of arguments can also be made for the second part of the Hamiltonian 
constraint T{N). The connection components Aa in (|169|) can be replaced by holonomies 
along edges with suitable tangents, and the integrated exterior curvature K by Poisson 
brackets of V with the regulated Euclidean part (I174p . as per (|17ip . The resulting ex- 
pression is quite complicated and contains ambiguities, but the correct refinement limit is 
obtained for a large class of regulators. 

Quantization 

We will now come to the quantization. The general idea is clear: Pick a family of regulators 
which converge to the continuum result. Replace Poisson brackets by commutators, and 
holonomies and volume operators by their operator counterparts, and obtain operators 

3 

cPiN) - ^E^(^n)E (('^MD) - ^-(D)) K\d) [hs,iO),v]). (175) 
n 7=1 

on the kinematic Hilbert space. Now take the refinement limit D — > E to obtain an 
operator Ce- There are, however, several difficulties when putting this program into 
practice: 

1. In the limit of infinite refinement, the operator is in danger of creating infinitely 
many loops and edges. Hence the limit may be ill defined. 

2. Even if problem 1. can be overcome, the operator will generically not converge, since 
typically 6*^°'* _L C^°'% for regulators D 7^ D'. 

3. Since h and V do not commute, there are ordering ambiguities. 

4. There is a lot of ambiguity in the choice of regulators since now there is no guaran- 
tee that different families of regulators will converge to the same operator, if they 
converge at all. 

The first problem can be solved by a suitable ordering. Let us consider the action on 
a spin network. The volume operator acts only at the vertices, hence ordering it to the 
right will force the loops and edges that are created by Ce to be attached to the vertices 
of the spin network only. Thus, for a given spin network, only finitely many new edges 
and loops can be created. This also partially solves problem 3. To deal with the rest of 
the difficulties, we will be less ambitious, and not demand convergence in the kinematic 
Hilbert space. Rather, we consider the matrix elements of Cg between one kinematic 
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state and one difFeomorphism invariant one. It turns out that due to the diffeomorphism 
invariance of the one state, many of the ambiguities in the attachment of the loops and 
edges do not change the matrix elements. What is more, for several types of regulators it 
is known that the matrix elements converge, 

lim (^IC^^Vt) is weU defined. (176) 

Typically, the matrix elements already become constant at a finite refinement, namely 
when the decomposition of E into cells is already so fine that there is at most one vertex 
of 7 per cell. 

Now we have to be careful. Convergence of the above matrix elements does not imply 
that there exists a limit operator on the kinematic Hilbert space. Rather, we can interpret 
(^|Cg I as an element in the (algebraic) dual space of Cyl, and hence conclude that there 
is an operator 

C'l : Hdiff -^ Cyl* . (177) 

The detailed features of this operator depend on the chosen family of regulators. But the 
generic features do not: 

• Cg acts locally at the vertices. 

• It acts by creating and annihilating edges and loops. 

One can proceed in the same way with the quantization of T(N), but since the quan- 
tized expression contains double commutators with Cg , the operator action becomes 
extremely complicated. Nevertheless it is well defined and finite. 

Solutions 

Given the definition of the Hamilton constraints we sketched above, what are the solutions? 
They are states ^ in "Hdiff such that 

(*|C(iV)/) = for all / G Cyl and all N. (178) 

One simple solution is the LQG vacuum | ), which can also be interpreted as a state in 
Hdiff- But more complicated solutions exist. For working out the set of solutions in some 
detail, details of the regularization used in the quantization of the constraints have to 
be fixed, since they do matter. Suffice it to say that so called exceptional edges play an 
important role in the construction of solutions. Exceptional edges are edges of the type 
created by the quantum constraint itself. We will not discuss this in detail, but refer to 
[731 [71I731 [S] for more detailed accounts. 

Solutions lie in the intersection of the kernels of all Hamilton constraints. Formally, 
the projector on this space can be expressed and approximated as follows 76 : 



Pc = 5{C) = I DiVe*^(^) 



l + i I DN I N{x)Cix) + ^-^ I ^N 1 1 N{xi)N{x2)Cixi)C{x2) + 



(179) 



C{x) denotes the local action of the constraint, which is zero unless x is the position of 
a vertex of the state acted upon. The path integral over N gives an infinite result, but 
by requiring diffromorphism invariance, it can be split into a divergent term that can be 
normalized away, and a finite remainder [77) . 

The matrix elements of the projector can then be expanded into a series 

oo 

{T,,\PcT,,) ^ J2 E'-'E^"— «(^7il^(«i)^(«2)...C(«^)|r^J (180) 



Ar=0 vi 
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where the finite sums are over all vertices of 72 and Cvi...vn is the finite remainder of the 
integral over the lapse function. It only depends on the diffeomorphism equivalence class 
of the vertex set {wi, W2, ■ • • wat}- We note that a priori the multiple applications of the 
local constraint in (jl80p do not make sense, since we have up to now only defined the 
constraint operators in such a way that domain and range are disjoint, see (|177|) . But it 
is possible to enlarge the domain of definition in such a way that multiple applications of 
the constraints become possible [771 [ZH] ■ We will sketch how this is done when we discuss 
the question of anomalies below. 

These matrix elements are interesting, because in principle they contain all the infor- 
mation about the inner product on the Hilbert space of physical states, 

{T,, \PcT,,) = {Pc T,, I Pc T^Jphys. (181) 

The expansion (|180p can be interpreted as a kind of Feynman expansion, organized in 
terms of how many times the constraint acts. The individual terms can be nonzero only 
if the action of the constraint operators on T^^ produces exactly T-^-^ . Thus the non- 
zero diagrams can be thought of as terms coming from the evolution of one spin network 
state into another. More precisely, they can be labeled by a two-complex, whose faces 
carry representations and whose edges carry intertwiners. The complex has the graphs 
71,72 as boundaries, and the internal vertices correspond to the cation of the constraints. 
These diagrams are called spin foams, and they show up independently in approaches that 
discretize the covariant path integral for gravity, compare the contribution by Rovelli. 
That they show up in an expression for the physical inner product of the canonical theory 
is a very encouraging link between canonical and covariant picture. In fact, in the light of 
the recent developments that are treated in Rovelli's contribution to this volume, we are 
getting close to actually having a precise correspondence 

quantum Hamilton constraint < — >■ spinfoam model. (182) 

We will now discuss some further aspects of the Hamilton constraint quantization. 

Symmetry, anomaly freeness, ambiguities 

In principle, it would be desirable to produce a symmetric, or even selfadjoint Hamiltonian 
constraint, 

C\N)^C{N). (183) 

But this turns out to be hard in practice, and there are even some no-go theorems [78] . 
Interestingly, there are heuristic arguments to the effect that one can not have both, 
symmetric constraints and a constraint algebra that is anomaly free. 

We have seen that the constraints classically close to form an algebra with respect to 
the Poisson bracket. The same should happen on the quantum level, now with respect to 
the commutators. Otherwise the gauge symmetries may have been broken when quantizing 
the theory. Such an anomaly in the gauge symmetries would strongly suggest the quantum 
theory to be unphysical. In particular, we are interested in the commutators 

[C{M),C{N)] (184) 

since by the above construction, we can already see that the Hamilton constraints trans- 
form correctly under gauge transformations and diffeomorphisms. Classically the above 
commutator is proportional to a diffeomorphism constraint, hence at minimum one re- 
quires that the commutator should vanish states of Hdiff- The problem is that the con- 
straints map Hdiff to a certain subspace of Cyl* which is strictly larger than T-Ldis- So the 
above commutator is not well defined, as it stands. There are two proposed solutions to 
this problem. The first, by Thiemann [75], is to look at the commutator on "Hkin, before 
removing the regulator. He finds 

[C^°\M),C^°\N)]^ something ^Q, (*| something = for |^) e Hdiff. (185) 
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In this sense, 

[C(M),C(iV)]|„,,, =0, (186) 

and the quantization is anomaly free. The other solution to defining the commutator is 
by Lewandowski and Marolf [7S] . They introduce a certain class of elements of Cyl* that 
is slightly larger than "Hdiff- Without going into technical details, a vertex- smooth state 
1^) is a state 

I*) e Cyl* : (*|C/0/^) is a function of V{(t){-f)), (187) 

i.e., of the set of vertices of the graph (jyil), for any diffeomorphism 0. Trivial examples of 
vertex smooth states are given by diffeomorphism invariant states. A less trivial example 
is the linear functional given by 

*' ^ (*| / iV^detgll*') (188) 

for a lapse function N and \^) in HdiS- 

Lewandowski and Marolf observe that {'^\C{N) is vertex-smooth for a large class of 
regulators, and that its action can be extended to vertex-smooth states. Moreover, they 
find 

(*vs|[C(M),C(7V)]-0, (189) 

where ^I'vs is vertex-smooth. As far as diffeomorphism invariant states are concerned, 
this result would be expected for an anomaly free representation. But since it holds for 
all vertex-smooth states, it is surprising and a little worrisome, since the term in the 
Dirac algebra that results from the Poisson bracket of two Hamiltonian constraints, a 
diffeomorphism constraint, would be expected to act non-trivial on most vertex smooth 
states. But this has to be checked explicitly, and it may be possible to find quantizations of 
this term that indeed vanish on vertex-smooth states. New light on this question may be 
shed by new results of Laddha and Varadarajan [TSl HOI US] , who employ new techniques 
to define constraints and their commutator algebra. 

We should not finish without pointing out that there are various ambiguities in the 
above procedure that are poorly understood, for example regarding the loop attachment 
and the representation of the newly created links (see however 181]). Overall, it is however 
very encouraging that we can find a family of well defined constraint operators that are 
anomaly free in a certain sense, and that lead to a convergence of the canonical and 
the spin-foam picture. Given the complexity of the Hamiltonian constraints of general 
relativity, these results are highly non-trivial. 

Some of the techniques discussed in this section will be used in the following two sections 
where the quantization of the two classical models introduced in section 12.2.11 and 12.2.21 

4.2 The quantum Brown-Kuchar model 

In this section we will discuss the quantization of the Brown-Kuchar model introduced in 
section 12.2.11 We want to quantize the reduced phase space whose elementary variables 
are given by the observables A^ and E-*, shown in equation (|29|) . Hence we need to look at 
the algebra of these elementary observables in order to know what kind of representations 
are possible for the corresponding quantum theory. A property of those models where 
deparametrization occurs is that the algebra of the elementary observables is isomorphic 
to the kinematical one, that is 

{A/(a),E^(a')} = 'lsiSpHa,a') (190) 

In general the algebra of observables can be more complicated and is given by the expres- 
sion [82] 

{OA,T,OE,T}^O^A,Er,T (191) 
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which involves the Dirac bracket denoted by {.,.}*. Here we denoted the general ob- 
servables associated with A, E with respect to a reference field T by Oa,t and Oe,t 
respectively. The reason why the Dirac bracket occurs on the righthand side is that the 
originally first class constraint of a given system together with the gauge-fixing constraints 
for the clock fields C/ := Tj — tj , where / labels the individual reference fields of a given 
model, form a system of second class constraints. Applied to general relativity, we would 
start with the first class system given by the Hamiltonian C and spatial diffeomorphism 
constraint Ca ■ Then for the Brown-Kuchar- model we choose four reference fields T, S^ 
with j = 1, 2, 3 and obtain four gauge-fixing constraints Co = T — t = 0, Cj = S^ — a^ = 0. 
In the framework of gauge unfixing introduced in [83] the construction of the observables 
in section 12.2.11 corresponds to transforming the second class constraints C and C'a again 
into first class constraints and using the gauge-fixing constraints to construct a projector 
that maps A and E onto their corresponding observables with respect to the now first class 
constraints C and Ca- Thus, the way how observables are constructed in the relational 
framework is a particular case of the gauge unfixing procedure. 

It seems that for the Brown-Kuchar-model the quantization of the reduced phase space 
seems to be a trivial task and considering only the algebra it looks like even a Fock quan- 
tization would be possible. However, this is not the case because likewise to the Dirac 
quantization where one requires that the kinematical representation needs to allow to 
implement the constraints as well defined operators, here we are only interested in those 
representations in which the physical Hamiltonian Hphys in p3|) can be implemented 
as a well defined operator. Since Hphys consists of terms that involve the gravitational 
contribution of the Hamiltonian and diffeomorphism constraint, Fock quantization is ex- 
cluded. However, a possible representation would be the one of the (gauge invariant) 
kinematical Hilbert space in LQG on which the constraint operators can be defined. Note 
that this representation becomes physically in this model since we are quantizing the re- 
duced phase space here. Hence, as a first possible representation for 'Hphys let us choose 
^phys = -^2 (-4, dfiAL), restricted to its gauge invariant subspace. Now, our task is to quan- 
tize the generator of the dynamics, that is Hphys. On the classical reduced phase space 
the expression h^(cr) = C^ — q^^CjC^ is constrained to be positive. Implementing this 
in the quantum theory would correspond to defining self-adjoint operators for h^(cr) and 
restrict for each a the spectral resolution of the Hilbert space to the positive part of the 
spectrum. Since this is technically impossible at the moment because of the complexity 
of the operators in the full theory, we use the absolute value under the square root and 
instead and quantize 

Hphys = /dV./|C2-q.'-fcCjCfc|(a) (192) 



5 

As a first step, likewise to the construction in section lTTl we need to regularize the classical 
expression. For this purpose we choose a partition of the spatial dust manifold S into 3 
dimensional cells D such that 5 = (J D. Hence, Hphys can be written as 

Hphys = E / rf''^v'|C2-qJ-^C,Cfc|(a) (193) 

Let us denote the refinement limit in which the partition becomes the continuum by 
D — >■ 5, the volume of the cells by Vo(n) and a point inside D by cr(n), then we can 
rewrite Hphys as a limit of a Riemann sum 



Hphys = ^im^^T/o(n)V|C2 - qJ-'=C,Cfe|(a(n)) (194) 

As a second step we will reformulate the expression under the square root above so that 
we are able to use quantization techniques that have been successfully applied to the case 
of the Hamiltonian constraint and have been discussed in section 14.11 For simplicity we 
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will restrict our discussion to the euclidean part of the Haniiltonian constraint Ce- As 
has been explained in section 14.11 once this part has been quantized the remaining part 
can be quantized using the operator for Ce- We introduce the (rescaled) magnetic field 
B-'j and its contraction with a co-triad given by 

BS := 2e^'=^F^„ B := B^rje, (195) 

where tj = —iajjl are, as before, a choice of a basis for su(2), gj are the Pauli matrices 
and F^j, is the curvature associated to A''. Note that the index position of the capital 
indices is not important here since these are the su(2) Lie algebra indices are pulled with 
(5''^. Using that TrijjTK) = —\^jk we obtain 

rr(B) = -eJ'"F)^^e/ (196) 

The co-triad can be expressed in terms of triads by the following formula 

2 det(e^^) 

together with the identities Ej = ■\/det(q)e-'^ and det(e'') = sgn(det (e) ) -y/det (q) we 
obtain 

^JKLt^J xpfc Tpi 

Tr(B) = -sgn(det(e)) J^ ^ (198) 

Vdet(q) 

Consequently we have [Tr(B)]^ = C'^. For the second term under the square root 
q^'^CjCfc using q^'*^ = E^jE^ ^^ / det{q) and C^ = FJ^E|- we obtain 

q^'^C.Cfc = ^ \ 199) 

^ ^ det(q) ^ ' 

On the other hand when we use that Tr{TjTjTK) — j^ijk and consider the term 4Tr(BTK) 
we obtain 

4Tr(BrK) = -sgn(det(e))^^M^^ (200) 

A/det(q) 



Thus we have 

Let us introduce the quantities 



q^^C.Cfc = 16Tr(Brj)Tr(BTK)5''^ -: -JJ^CjCk (201) 



C(n) := j d?aC{a) Cj{U) - f d^aCj{a) 



a D 

then in the refinement limit we can rewrite Hphys as 



(202) 



Hphys - ^im V Vic^(n) - s-'KCj{n)CK{a)\ (203) 

and this finishes the regularization of the classical expression. Using the notation f^ := 
(— I2,4tj) with /Lt = 0, 1,2,3 and the classical identity in (|167p we can rewrite the regu- 
larized expressions as 

C^{U) = j d?aTr{BT^) = ^ j Tr{¥ A {W{U),A}r^) (204) 

D D 

where 1^(0) is the volume of D given by 

y(n) = /dVv/det(q) (205) 
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The corresponding quantum operator is then defined as 



Hphys = hm V W cS(n)Co(n) - sJ'<&jiD)CKia) (206) 

a^s ^ V I 

and one needs to show that the hmit yields a weU defined operator on "Hphys- For the 
operator Co(n) this has been shown in [26| and is briefiy discussed in section l4Jl The 
operator C/(n) can be quantized using similar techniques since in its definition also enters 
the Tr(B) term with an additional tj matrix inside the trace. At this point the symmetries 
of the classical physical Hamiltonian become important. As mentioned in section [2.2.11 
the physical Hamiltonian is invariant under (active) diffeomorphisms on the dust manifold 
S and we would like to preserve this symmetry also in the quantum theory having the 
consequence that Hphys needs to be implemented as a spatially diffeomorphism invariant 
operator. The representation we choose for Hphys is the gauge invariant sector of the usual 
kinematical representation of LQG, namely L2{A,d^AL)- For this representation it was 
shown in 44 that spatially diffeomorphism invariant operators, need to be quantized in a 
graph preserving way, meaning that those operators do not modify the underlying graph 
that a spin network function is defined on. However, in its usual quantization discussed 
in section 14.11 the operator Co is quantized in a graph changing way. Thus, choosing 
the usual kinematical representation of LQG for Hphys together with the requirement 
that the classical symmetries of Hphys carry over to the quantum theory forces us to 
quantize the operators C^ in a graph preserving way. Since the Hilbert space in the chosen 
representation decomposes into an orthogonal sum of the Hilbert spaces associated with 
each individual graph, this means that each of these graph Hilbert spaces needs to be 
preserved separately. In order to implement this graph preserving property we introduce 
the notion of a minimal loop: Given a graph 7, consider a vertex v € 1^(7) and a pair of 
edges e, e e ^^(7) of edges starting at the vertex v. A loop a^.e,g in 7 starting at v going 
a long the edge e and ending at v along the edge e~^ is said to be minimal provided that 
there exist no other loop in 7 with these properties and fewer edges transversed. Using 
the notion of a minimal loop we can define an operator for each graph at a given vertex v 

>\-^-^^>\ (ei,e„e3)GT„(7) l^7.-.e..e.| „gi^_^^^_^^ 

(207) 
here Ty{'-f) denotes the set of ordered triples of edges at the vertex v, whose tangent vectors 
at V are linearly independent, and L^yf,i,ej is the set of minimal loopf|j Furthermore, 
Vry^v is the Ashtekar-Lewandowski volume operator shown in equation (j88p . The operator 
for the physical Hamiltonian for each graph 7 is then defined as 



Hphys,7= Y. 'J\Pi\p^^-^i,--^"'^ln,v^K,^,v)P^\ (208) 

i>GV(7) 

whereby P-y : Hphys — > ^phys 7 ^^ ^^ orthogonal projection operator defined analogously 
to the orthogonal decomposition shown in (j69p . that needs to be introduced in order to 
ensure that Hphys, 7 is graph preserving. Although the loop is attached along already 
existing edges of the graph it can still be the case that when a holonomy operator is 
acting the resulting product of representation includes the trivial one. The final operator 
is then defined as 

-Hphys ~- \\j phySi7 \Z\j\)) 



Note that the introduction of the minimal loop here causes no further complication when the semiclassical 
limit of this operator is considered. This can for instance be seen in jSJ] where the semiclassical limit in the 
context of AQG is discussed and also the concept of the minimal loop is used in the quantization. 
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The fact that Hphys has to be quantized in a graph preserving way conies from the spatial 
diffcomorphism invariant corresponding classical expression of the physical Hamiltonian. 
However, one might take the point of view that including those projection operators into 
the operator Hphys to enforce the graph preserving property of the operator looks slightly 
artificial. A way to avoid this issue and thus also the projection operators in (|208p is 
to change the representation for Hphys- One possible other representation introduced 
in the framework of Algebraic Quantum Gravity (AQG) |29| is von Neumann's infinite 
tensor product representation (ITP). In the context of AQG one does not work with the 
embedded graphs used in LQG but considers one (fundamental) abstract combinatorial 
graph on which the quantum dynamics is defined. The embedding of the graph into a 
given spatial manifold happens only in the semiclassical sector of the theory and how 
the abstract graph is embedded is encoded in semiclassical states. To each edge of the 
abstract graph one associates an L2{SU{2),d^H) Hilbert space and one considers a graph 
with countable infinitely many edges. One of the motivations to introduce the AQG model 
was that semiclassical computations of dynamical operators technically simplify in this 
setup. This is due to the fact that for graph preserving operators the current existing 
semiclassical states can be used and those operators can be defined more naturally in 
the AQG framework. In [M] a combinatorial graph of cubic topology was chosen and 
considering this graph we can define an (algebraic) operator for the physical Hamiltonian. 
The algebraic version of the operator in (|207l) is given by 

C^,.= E *i*2S3e^^^^'^Tr(f^/i„,^^,^,,^^^/ie.,.3.3fc4^^,K]) (210) 

here e^js denotes the edge starting at v and going in positive (s=+l) or negative (s=-l) 
I-direction and ajsjs is the unique minimal loop in the algebraic graph of cubic topology, 
that starts at v goes along the edge ejg and comes along the edge ej^ before ending in v 
again. In analogy one can define an algebraic volume operator, that is for a cubic graph 
of the form 



V — /^ 



^ 2^ e eLMNJe^Jsi'Je^,Js2'^e^,Ks3 

Sl,S2,S3 = ±l 



(211) 



The final algebraic operator is then given by 

(212) 



-Hphys — 7 ^ 

vev{j} 



ClC,-S-'Kc\C 



J,v'^K,-, 



where the sum runs over the countable infinitely many vertices of the algebraic graph. In 
contrast to the LQG framework in AQG trivial representations associated to the edges 
are allowed. The picture of the dynamics is then that dynamical operators do not change 
the underlying infinite abstract algebraic graph but only representations associated to the 
edges. However, since trivial representations are allowed subgraphs of the fundamental 
algebraic graph can and will be modified so that the quantum dynamics in the algebraic 
framework looks similar to the graph modifying one in LQG. In both formulations the 
usual LQG and the AQG model geometrical operators are defined on the physical Hilbert 
space and are thus observables. Hence, these models are two examples where the discrete 
spectra of those operators is carried over to the physical sector. 

A further possibility to formulate a model for LQG, in which operators can be defined in a 
graph changing way is the scalar field model discussed in section r2.2.2l whose quantization 
will be discussed in the next section. 

4.3 The quantum scalar-field model 

In this section we discuss the quantization of the scalar field model whose classical theory 
was introduced in section 12.2.21 and we will closely follow the presentation from [SF . As 
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already mentioned at the end of section 12.2.21 the physical Hilbert space "Hphys of this 
model will be constructed from the gauge invariant subspace of the diffeomorphism in- 
variant Hilbert space denoted by H^-^g for the reason that the diffeomorphism as well as 
the Gauss constraint are solved by means of Dirac quantization. As discussed in section 
13.3.21 and 13.3.31 the diffeomorphism invariant Hilbert space can be constructed by using a 
rigging map rj^ig : 2?kin -^ ^diff ^ -^kin- ^^ addition in order to construct the operator 
corresponding to the physical Hamiltonian i/phys of this model, we need another Hilbert 
space denoted by 'Hdiff,x associated to a subgroup denoted by Diff{T,,x) of Diff{'E) 
including those diffeomorphisms, which preserve a given point a; G S. The construction 
of the rigging map for Di//(E, x) works analogously to that of Dif f{Yj) when using the 
techniques introduced in section 13.3.21 Next the gauge invariant subspace of these two 
Hilbert spaces can be easily obtained because for each gauge invariant cylindrical function 
/ G Cyl, the linear functional 77diff(/) is not affected by gauge transformations that act 
on / G Cyl in the sense that the expression [?7diff (/)](/) is invariant under gauge transfor- 
mation. Consequently, we obtain "H^jff and Ti^g ^ respectively by restricting I?kin =Cyl 
to the subspace of gauge invariant cylindrical functions. 

When we decide to do not reduce with respect to the Hamiltonian constraint at the clas- 
sical level, we need to construct solutions to the Hamiltonian constraint in the quantum 
theory, this yields an equation of the form 

(tt{x) - h{x)'\ * = (213) 

Taking into account that tt is quantized as --iS/6(j){x) (setting ?i = 1) we obtain as a 
(formal) general solution 

*((/), A) = eJ'^ d^^4>ix)hix)^^^-^ (214) 

where '0(A) is an SU(2) gauge and spatially diffeomorphism invariant function. 

Physical operators, these are operators that correspond to classical Dirac observables, 
wiU be defined on "Hphys or a dense subspace X'phys of it. We explained in section [3.3.21 
that the action of operators on T^kin can be extended to T)^^^. Therefore, in our case, we 
have for symmetric diffeomorphism and gauge invariant operators L defined originally on 
"Hkin a natural action on Hdiff, where we denote the extended operator by L' , given by 

[^''7difF (/)](/) := [VdmimL^f) = [^diff(/)](i/) = (^diff(i/) , ^diff(/))difF (215) 

So far, the operators L are only observables with respect to the diffeomorphism and 
Gauss constraint but not with respect to the Hamiltonian constraint. In this model those 
observables are not constructed at the classical level, like in the Brown-Kuchar model of 
section 12.2.11 but are constructed as Dirac observables directly in the quantum theory on 
"Hphys- A quantum Dirac observable O is defined as an operator on 'Hphys (or a dense 
subspace X'phys) with the following properties 

• O is SU(2)-gauge and spatially diffeomorphism invariant. 

• The operator O commutes with the Hamiltonian constraints, that is [C^°*'{x), O] = 
for all x G E. 

Inspired by the relational framework [30l [31] for the classical theory, one can (formally) de- 
fine a family of Dirac observables. Let L be an SU(2) gauge and diffeomorphism invariant 
linear operators, then the operator 0{L) defined as 

i fd x( d>(x) — T(x))h(x) ^ —i fd x( d)(x)—T(x)]h(x) 

is a quantum Dirac observable. Here r is the value that the reference field (p takes while 
being transformed along its gauge orbit. We need to choose r(x) :— t with r being 
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a constant real number in order to ensure that the resulting operator is spatially dif- 
fcomorphism invariant. Note that the expression above can also be obtained from the 
requirement that the (formal) solutions shown in (|214p should be mapped into solutions 
by Dirac observables. The classical interpretation of these Dirac observables is precisely 
the analogue of those formal power series used to construct the observables A/ and E-* 
in the Brown-Kuchar model shown in equation 
theory 



.J J 

Explicitly, we have in the quantum 



Or{L)=Y.'- 

'—^ Til 



i, Wx)-T)h{x)) 



(n) 



which is the quantization of the classical expression 



ra=0 ■ ■l 



(217) 



(218) 



The quantum dynamics of the observables is given by the following equation 

d 



l — Or{L) 

dr 



0(i),-ffphys 



(219) 



and those equations are the analogue of the Heisenberg picture in quantum mechanics. 
Our final task is to implement the operator -ffphys that is a quantization of the classical 
expression in p9|) . Promoting the individual terms under the square root in Hphys to 
operators we obtain the heuristic expression 



H, 



phys 



= d 



C2 - q-bCaCb 



(220) 



where q := det(g). This operator will be defined on a suitable domain of H^jfj and thus 
act only on states that are spatially diffeomorphism invariant as long as we restrict to 
the operator evolution defined by ()217|) . The operators corresponding to the classical ex- 
pression q°'^CaCb should annihilate diffeomorphism invariant states. Therefore, assuming 
a suitable operator ordering for Hphys we assume that we can work with the simplified 
operator 

d^xh{x) (221) 



H, 



phys 



d^a 



and this is also the physical Hamiltonian suggested by Rovelli and Smolin in [3S] . Although 
we can use some of the already existing quantization techniques in the literature for 
C and \/q respectively, what we need to define is an operator for y^~2y/qC and the 
already existing operators are the Hamiltonian constraint C smeared against arbitrary 
lapse functions (graph- modifying) in [26 , the master constraint mentioned in section [2.21 
[2 7) (graph-modifying and graph-preserving) and the physical Hamiltonian of the Brown- 
Kuchar model (graph preserving) in j34| . The physical Hamiltonian i?phys for the scalar 
field model will be defined on a suitable domain of "H^jg in a graph changing way. Let us 
briefly sketch the quantization procedure. The regularization of the classical expression 
for C will be that of [5], where the operator valued distribution for C is defined as 



d^xN{x)C{x) 



xen 



Nix)C', 



(222) 



One of the differences to the regularization chosen by Thiemann in |26j is the way how 
the loop is attached to the graph on which the spin network functions are defined on. In 
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Thiemann's proposal the loop runs along the edges e/ and ej that belong to the graph 7 
(see section 14. ip whereas in the regularization in [S] the loop lies in the plane spanned by 
the edges e/ and ej but is only connected to the graph at the vertex v. Note that neither 
the operator ^C{x) nor the operator h{x) are defined on ^Jfig due to their dependence 
on X that breaks difFeomorphism invariance. Therefore, we need the Hilbert spaces "H^jg ^ 
in order to implement those operators in the quantum theory. Each individual operator 
C'^ maps its domain from Hdiff to Hdiff.x and as discussed in [5S] defines naturally an 
operator on 'Hdiff.x- The operators C'^ from [5] are non-symmetric, which is no problem 
as long as we are working with constraint operators. However, here we are quantizing a 
physical Hamiltonian, that involves a square root and we will need to perform a spectral 
decomposition in order project on the positive part of the spectrum of the operator under 
the square root. Therefore, we consider the operators 

C.:=1{CL+C'J) (223) 

The operator valued distribution, that we consider for implementing the Hamiltonian 
density h{x) of -ffphys is then of the form 

C{x) := Y^ S{x,x')C^, (224) 

x'en 

Likewise one can define an operator valued distribution for the classical expression ^/q 
given by 

Vt(a:)= ^5(x,x')yt,, (225) 

Note that both operator values distributions above are well defined for the reason that 
when they are smeared against an arbitrary smearing function F and applied to some 
cylindrical function only a finite number of non zero terms occur in the sum 

/• _ ^ _ 

/ d^xF{x)^{x)f^ = Y. PMVq.J-f (226) 

E fc=l 

where vi,...,vn are the vertices of the graph 7 and similar for the operator C{x). The 
expressions in (j224|) and (j225|) can then be used to define an operator for ifphys — /^ h{x) 
on (a subspace of) T-L^^g ^ 



h{x) := Y. Six,x')^-2^1C,,^1 (227) 

Due to the square root in the equation above h{x) is only well defined on the subspace 

of H^iff X where the spectrum of ^ ,Cx'y/q', is positive. In order to be able to consider 
only the positive part of the spectrum we need to choose a selfadjoint extension for h{x) 
and in general this choice might not be unique. Let us denote the subspace of Hdiff.x 
corresponding to the positive part of the spectrum by Hdiff,x,+- As discussed in [85] there 
exists a natural map 77s : "Hjfiff ^ -> "H^jfj with 77diff(E,x)(/) ^ '7diff(s)(/)- The domain of 
the physical Hamiltonian operator is then the image of 77s on Hdiff,x,+ and we obtain as 
the physical Hilbert space Hphys = ??E(^difr,x,+) with inner product 

(e/. d'^.fh^^ g/, d^-4>h^')^^^^ := (^, ^')^,g (228) 

with t/j,ip' G "H^iff- The final form of the operator for iJphys is then given by 



Hphys = / k^) = Y V -2vt!6. vt! (229) 
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Finally, let us comment on the classical symmetries of the Hamiltonian in this m.odel. 
Likewise to the Brown-Kuchar model iJphys is invariant under spatial diffeomorphisms. 
However, for the reason that here a different representation than the kinematical rep- 
resentation of LQG was chosen, the requirement that iJphys needs to be quantized in 
a graph preserving way is absent. Furthermore, also the classical Hamiltonian densi- 
ties {h{x),h{y)} = commute a general property of those deparametrized models. In 
the framework of the habitat, being the home of so called vertex smooth states, briefly 
mentioned at the end of section 14. 1[ it was shown in [7S] that the commutator of two 
Hamiltonian constraints smeared against arbitrary lapse functions vanishes on the habi- 
tat. Thus, one would expect that the Hamiltonian densities in the scalar field model will 
also vanish on the habitat [h{x), h{y)] = 0. 

Although the graph modifying property of the physical Hamiltonian in the scalar field 
model has some advantages as far as the rather artificial infinitely many conservation 
laws for graph preserving operators in the Brown-Kuchar model are concerned, when the 
semiclassical sector of the model is considered the analysis becomes more complicated. 
The reason for this is that with the current existing semiclassical techniques for the full 
theory, graph modifying operators cannot be analyzed and those techniques would need 
a strong improvement in order to be able to analyze the semiclassical sector of the scalar 
field model in full generality. 

Summarizing, both quantum models, the Brown-Kuchar as well as the scalar field model 
could be a step in the right direction and help in the future to extract some physical 
information out of full LQG but of course this is a long range project and on the way 
there are still many technical issues to solve. 

5 Summary and open problems 

In the present lectures we have given an introduction to loop quantum gravity, in particular 
its canonical quantization techniques. We also saw how it makes contact with the path 
integral formulation developed in spin foam gravity. The latter topic, and its connection to 
what we presented here is covered in detail in the contribution by Rovelli. But also in what 
we did cover, we have left out many details, and did not even touch on any applications, 
such as to the quantum theory of black hole horizons (covered in the contribution by 
Barbero, Lewandowski and Villasenor), or to cosmology (covered in the contribution by 
Singh). A good starting point to get an overview over all of these developments is the 
contribution by Ashtekar in these proceedings. 

However, we hope that we have explained at least some of the other big achievements of 
loop quantum gravity, namely its description of quantum geometry and the corresponding 
dynamics. The quantum theory of (extrinsic and intrinsic) geometry, as described in 
section[2]comprises in particular geometric operators with a discrete spectrum, the scale of 
which is set by Planck lenght, and diffeomorphism invariant states. We have furthermore 
seen that, based on this, well defined Hamiltonian constraints, and in the case of the 
matter models we considered, well defined Hamiltonian operators can be obtained. This 
is a highly non-trivial result, given the complicated nature of the classical dynamics. 
Moreover, there is a clear connection to the spin foam approach to loop quantum gravity. 

Although many of the structures that we have described have already been investigated 
for some time, there are still lots of new developments. Some of those have already 
been mentioned in the main text, but there are many more that we could not cover 
in these lectures, among them new connections between the full theory and symmetry 
reduced models [86l [87] , coherent states for quantum geometry [88l HH |90l EU |92l [93] , an 
interpretation of quantum geometry in terms of polyhedra in flat space [Ml [Ml [Ml [Ml 
[97[ [98] , and corresponding stunning results about the quantum volume [99] and the use 
of spinor techniques iMl \IM [Ml ^M [103] . 

Let us finally list some important questions that are the subject of ongoing investigation 
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in loop quantum gravity: 

• Barbero-Inimirzi parameter: What role does it ultimately play in loop quantum 
gravity with and without matter? 

• Controlled approximations: Loop quantum gravity is a non-perturbative approach 
to the quantization of gravity, but approximations will be vital to do physics. How 
can we find controlled approximations to situations with symmetries from the full 
theory, or approximately solve the Hamilton constraints? In the context of the 
reduced models the Hamiltonian constraint is already solved and what we end up is 
an evolution equation in the physical Hilbert space. However, also here one needs 
approximation techniques for the reason that the evolution equations are similar 
complicated to the solution equations of the Hamiltonian constraint. 

• Loop quantum gravity and matter: Which types of matter can be consistently cou- 
pled to loop quantum gravity? What are the implications of quantized space-time 
geometry to the propagation of matter? 

• Physics from Hamilton constraints and Hamiltonians: How does one extract physics 
from the solutions to the constraints? In particular one should be able to understand 
how ordinary quantum field theory and classical general relativity are embedded into 
loop quantum gravity. The first should correspond to a sector of quantum gravity 
where quantum fluctuations of the geometry are small but matter is still treated 
as a quantum object, whereas for general relativity both the matter and geometry 
quantum fluctuations are expected to be negligible. Furthermore, it is important to 
analyze how ambiguities in the quantization of constraints and physical Hamiltonians 
do reflect in physical properties of the theory. 

• Connection to spin foam gravity: What is the precise relation between scattering 
amplitudes and the physical inner product? Which quantization of the Hamilton 
constraint corresponds to which vertex amplitude? 

For some of these, there are already some insights. Answers to these questions will be 
crucial for the path that loop quantum gravity takes in the future. 

Acknowledgments 

The authors would like to thank the organizers of the 3rd quantum geometry and quantum 
gravity school in Zakopane for their kind invitation. The authors would also like to thank 
Jerzy Lewandowski for carefully reading our review and giving valuable feedback that 
improved the current version of the review. This work was supported by the European 
Science Foundation through its network Quantum Geometry and Quantum Gravity. 

References 

[1] C. Rovelli, Quantum gravity, Cambridge, UK: Cambridge Univ. Pr. (2004) 455 p 

[2] T. Thiemann, Modern canonical quantum general relativity, Cambridge, UK: Cam- 
bridge Umv. Pr. (2007) 819 p |gr-qc/01100"3^ . 

[3] R. Gambini, J. PuUin A first course in Loop Quantum Gravity Oxford Univ. Pr. 
(2011) 183p 

[4] M. Bojowald, Canonical Gravity and Applications: Cosmology, Black Holes, and 
Quantum Gravity Cambridge, UK: Cambridge Univ. Pr. (2010) 312p 

[5] A. Ashtekar and J. Lewandowski, "Background independent quantum gravity: A 
Status report," Class. Quant. Grav. 21, R53 (2004) |arXiv:g r-qc/0404018, . 

[6] A. Perez, "Loop quantum gravity," Europhys. News 37N3, 17 (2006). 



50 



[7] A. Perez, "Loop quantum gravity: An introduction," AIP Conf. Proc. 1132, 386 
(2009). 

[8] P. Dona and S. Speziale, "Introductory lectures to loop quantum gravity," 
arXi v: 1007. 0402 [gr-qc]. 

[9] H. Sahlmann, "Loop Quantum Gravity - A Short Review," larXiv:1001. 41881 [gr-qc]. 

[10] A. Ashtekar, "New Variables for Classical and Quantum Gravity," Phys. Rev. Lett. 
57, 2244 (1986). 

[11] J. F. Barbero G., "Real Ashtekar variables for Lorentzian signature space times," 
Phys. Rev. D 51, 5507 (1995) arXiv:gr-qc/94100l4 . 

[12] G. Immirzi, "Real and complex connections for canonical gravity," Class. Quant. 
Grav. 14, L177 (1997) [arXiv:gr-qc /9612030|. 

[13] S. Hoist, "Barbero's Hamiltonian derived from a generalized Hilbert-Palatini ac- 
tion," Phys. Rev. D 53, 5966 (1996) [arXiv:gr-qc/9511026 . 

[14] A. Perez and C. Rovelli, "Physical effects of the Immirzi parameter," Phys. Rev. D 
73, 044013 (2006) jarXiv:gr-qc/0505081| . 

[15] S. Carlip, Quantum gravity in 2+1 dimensions, Cambridge, UK: Univ. Pr. (1998) 
276 p 

[16] N. Bodendorfer, T. Thiemann and A. Thurn, "New Variables for Classical and 
Quantum Gravity in all Dimensions I. Hamiltonian Analysis," larXiv: 1105.37031 [gr- 
qc]. 

N. Bodendorfer, T. Thiemann and A. Thurn, "New Variables for Classical and 
Quantum Gravity in all Dimensions II. Lagrangian Analysis," larXiv: 1105.37041 [gr- 
qc]. 

N. Bodendorfer, T. Thiemann and A. Thurn, "New Variables for Classical and 
Quantum Gravity in all Dimensions HI. Quantum Theory," larXiv:1105 .3705 [gr- 
qc]. 

N. Bodendorfer, T. Thiemann and A. Thurn, "New Variables for Classical and 
Quantum Gravity in all Dimensions IV. Matter Coupling," arXiv:1105.3706 [gr-qc]. 

[17] S. Alexandrov, "S0(4,C) covariant Ashtekar-Barbero gravity and the Immirzi pa- 
rameter," Class. Quant. Grav. 17, 4255 (2000) ! arXiv:gr-qc/0005085| . 

[18] S. Alexandrov, "On choice of connection in loop quantum gravity," Phys. Rev. D 
65, 024011 (2002) |arXiv:gr-qc/010707l. 

[19] S. Alexandrov and E. R. Livine, "SU(2) loop quantum gravity seen from covariant 
theory," Phys. Rev. D 67, 044009 (2003) |arXiv:gr-qc/0209105| . 

[20] E. R. Livine, "Projected spin networks for Lorentz connection: Linking spin foams 
and loop gravity," Class. Quant. Grav. 19, 5525 (2002) ,,arXiv:gr-qc/02070 84 . 

[21] F. Cianfrani and G. Montani, "Towards Loop Quantum Gravity without the time 
gauge," Phys. Rev. Lett. 102, 091301 (2009) jarXiv:0811. 19161 [gr-qc]]. 

[22] F. Cianfrani and G. Montani, "The Role of Time Gauge in Quantizing Gravity," 
■arXiv:0 904.0573 [gr-qc]. 

[23] F. Cianfrani and G. Montani, "The Immirzi parameter from an external scalar field," 
Phys. Rev. D 80, 084040 (2009) |arXiv:0907.1530l [gr-qc]]. 

[24] M. Geiller, M. Lachieze-Rey, K. Noui and F. Sardelh, "A Lorentz-Covariant Con- 
nection for Canonical Gravity," SIGMA 7, 083 (2011) ^rXiv: 1103.4057 [gr-qc]]. 

[25] M. Geiller, M. Lachieze-Rey and K. Noui, "A new look at Lorentz-Covariant Loop 
Quantum Gravity," Phys. Rev. D 84, 044002 (2011) |arXiv:1105.4194l [gr-qc]]. 

[26] T. Thiemann, "Anomaly - free formulation of nonperturbative, four-dimensional 
Lorentzian quantum gravity," Phys. Lett. B 380, 257 (1996) ' arXiv:gr-qc/9606088| . 

51 



[27] T. Thiemann, "The Phoenix project: Master constraint program for loop quantum 
gravity," Class. Quant. Grav. 23, 2211 (2006) |arXiv:gr-q c/0305080] . 

[28] B. Dittrich and T. Thiemann, "Testing the master constraint programme for loop 
quantum gravity. I. General framework," Class. Quant. Grav. 23, 1025 (2006) 
[arXiv:gr-qc/0411138| . 

B. Dittrich and T. Thiemann, "Testing the master constraint programme for loop 
quantum gravity. II. Finite dimensional systems," Class. Quant. Grav. 23, 1067 
(2006) [arXiv:gr-qc/0411139| . 

B. Dittrich and T. Thiemann, "Testing the master constraint programme for 
loop quantum gravity. III. SL(2,R) models," Class. Quant. Grav. 23, 1089 (2006) 
[arXiv:gr-qc/0411140. 

B. Dittrich and T. Thiemann, "Testing the master constraint programme for loop 
quantum gravity. IV. Free field theories," Class. Quant. Grav. 23, 1121 (2006) 
]arXiv:gr-qc/041114l . 

B. Dittrich and T. Thiemann, "Testing the master constraint programme for 
loop quantum gravity. V. field theories," Class. Quant. Grav. 23, 1143 (2006) 

arXiv:gr-qc/0411142 . 

[29] K. Giesel and T. Thiemann, "Algebraic Quantum Gravity (AQG). I. Conceptual 
Setup," Class. Quant. Grav. 24, 2465 (2007) |arXiv:gr-q c/0607099 . 

[30] C. Rovelh, "QUANTUM REFERENCE SYSTEMS," Class. Quant. Grav. 8, 317 
(1991). 

C. Rovelh, "TIME IN QUANTUM GRAVITY: PHYSICS BEYOND THE 
SCHRODINGER REGIME," Phys. Rev. D 43, 442 (1991). 

C. RoveUi, "QUANTUM MECHANICS WITHOUT TIME: A MODEL," Phys. Rev. 
D 42, 2638 (1990). 

[31] B. Dittrich, "Partial and complete observables for Hamiltonian constrained sys- 
tems," Gen. Rel. Grav. 39, 1891 (2007) 'arXiv:gr-qc/0411013". 
B. Dittrich, "Partial and complete observables for canonical general relativity," 
Class. Quant. Grav. 23, 6155 (2006) arXiv:gr-qc/ 0507106 . 

[32] J. D. Brown and K. V. Kuchar, "Dust as a standard of space and time in canonical 
quantum gravity," Phys. Rev. D 51, 5600 (1995) [;arXiv:gr-qc/940900l . 

[33] K. Giesel, S. Hofmann, T. Thiemann and O. Winkler, "Manifestly Gauge-Invariant 
General Relativistic Perturbation Theory. I. Foundations," Class. Quant. Grav. 27, 
055005 (2010) |'arXiv:0711.0115l [gr-ac]]. 

[34] K. Giesel and T. Thiemann, "Algebraic quantum gravity (AQG). IV. Reduced phase 
space quantisation of loop quantum gravity," Class. Quant. Grav. 27, 175009 (2010) 
[arXiv:0711.0TT9 [gr-qc]]. 

[35] C. Rovelh and L. Smolin, "The Physical Hamiltonian in nonperturbative quantum 
gravity," Phys. Rev. Lett. 72, 446 (1994) pXiv:gr-qc/9308002|. 

[36] K. V. Kuchar and J. D. Romano, "Gravitational constraints which generate a lie 
algebra," Phys. Rev. D 51, 5579 (1995) |arXiv:gr-qc/"95 6l005 . 

[37] A. Ashtekar, T. Pawlowski and P. Singh, "Quantum nature of the big bang," Phys. 
Rev. Lett. 96, 141301 (2006) 'arXiv:gr-qc/0602086|. 

A. Ashtekar, T. Pawlowski and P. Singh, "Quantum Nature of the Big Bang: 
An Analytical and N umerical Investigation. I," Phys. Rev. D 73, 124038 (2006) 
'arXiv:gr-qc/0604013]. 

A. Ashtekar, T. Pawlowski and P. Singh, "Quantum Nature of the Big Bang: Im- 
proved dynamics," Phys. Rev. D 74, 084003 (2006) l arXiv:gr-qc /0607039|. 

[38] V. Husain and T. Pawlowski, "Dust reference frame in quantum cosmology," Class. 
Quant. Grav. 28, 225014 (2011) a rXi^:1108.1147i [gr-qc]]. 

52 



[39 
[40 

[41 
[42 
[43 
[44 

[45 

[46 

[47 

[48 
[49; 
[50 

[51 

[52; 

[53 



C. Rovelli and L. Smolin, "Loop Space Representation of Quantum General Rela- 
tivity," Nucl. Phys. B 331, 80 (1990). 

A. Ashtekar and C. J. Isham, "Inequivalent observable algebras: Another ambiguity 
in field quantization," Phys. Lett. B 274, 393 (1992). 

A. Ashtekar and J. Lewandowski, "Representation theory of analytic holonomy C* 



algebras," gr-qc/9311010 



A. Ashtekar and J. Lewandowski, "Projective techniques and functional integration 
for gauge theories," J. Math. Phys. 36 (1995) 2170 gr-qc/9411G46, . 

A. Ashtekar and J. Lewandowski, "Differential geometry on the space of connections 
via graphs and projective hmits," J. Geom. Phys. 17 (1995) 191 [ hep-th/9412073| . 

A. Ashtekar, J. Lewandowski, D. Marolf, J. Mourao and T. Thiemann, "Quan- 
tization of diffeomorphism invariant theories of connections with local degrees of 
freedom," J. Math. Phys. 36, 6456 (1995) |arXiv:gr-qc/9504018] . 

J. Lewandowski, A. Okolow, H. Sahlmann and T. Thiemann, "Uniqueness of diffeo- 
morphism invariant states on holonomy-flux algebras," Commun. Math. Phys. 267, 
703 (2006) [arXiv:gr -qc/0504147. 

A. Ashtekar and J. Lewandowski, "Quantum theory of geometry. 1: Area operators," 
Class. Quant. Grav. 14, A55 (1997) |arXiv:gr-qc/9602046| . 

H. Sahlmann, "Some results concerning the representation theory of the al- 
gebra underlying loop quantum gravity," J. Math. Phys. 52, 012502 (2011) 
[arXiv:gr-qc/0207lTi] . 

A. Ashtekar, J. Lewandowski and H. Sahlmann, "Polymer and Fock representations 
for a scalar field," Class. Quant. Grav. 20, LU (2003) jarXiv:gr-qc/0211012j . 

L. Smolin, "Recent developments in nonperturbative quantum gravity," 
|arXiv:hep-t h/9202022, 

C. Rovelli and L. Smolin, "Discreteness of area and volume in quantum 
gravity," Nucl. Phys. B 442, 593 (1995) [Erratum-ibid. B 456, 753 (1995)] 
arXiv:gr-qc/9411005 . 

A. Ashtekar and J. Lewandowski, "Quantum theory of geometry. 2. Volume opera- 
tors," Adv. Theor. Math. Phys. 1, 388 (1998) [arXiv:gr-qc/9711031| . 

J. Brunnemann and T. Thiemann, "Simplification of the spectral analysis of the 
volume operator in loop quantum gravity," Class. Quant. Grav. 23, 1289 (2006) 
[arXiv:gr-qc/04050"6^ . 

J. Brunnemann and D. Rideout, "Properties of the volume operator in loop quantum 
gravity. I. Resuhs," Class. Quant. Grav. 25, 065001 (2008) arXiv:0706.0469 [gr-qc]]. 
J. Brunneman and D. Rideout, "Properties of the volume operator in loop quan- 
tum gravity. IL Detailed presentation," Class. Quant. Grav. 25, 065002 (2008) 
>rXiv:0706.0382' [gr-qc]]. 

J. Brunnemann and D. Rideout, "Oriented Matroids - Combinatorial Structures 
Underlying Loop Quantum Gravity," Class. Quant. Grav. 27, 205008 (2010) 
|arXiv:1003.2348l [gr-qc]]. 

[54] K. Giesel and T. Thiemann, "Consistency check on volume and triad operator 
quantisation in loop quantum gravity. I," Class. Quant. Grav. 23, 5667 (2006) 
|arXiv:gr-qc/050703^ . 

K. Giesel and T. Thiemann, "Consistency check on volume and triad operator 
quantisation in loop quantum gravity. II," Class. Quant. Grav. 23, 5693 (2006) 
[arXiv:gr-qc/"050 7037| . 

[55] T. Thiemann, "A Length operator for canonical quantum gravity," J. Math. Phys. 
39, 3372 (1998) arXiv:gr-qc/9606092| . 

53 



[56] E. Bianchi, "The Length operator in Loop Quantum Gravity," Nucl. Phys. B 807, 
591 (2009) |arXiv:0806.4710l [gr-ac]]. 

[57] Y. Ma, C. Soo and J. Yang, "New length operator for loop quantum gravity," Phys. 
Rev. D 81, 124026 (2010) [arX iv: 1004. 1063 [gr-qc]]. 

[58] D. Giulini and D. Marolf, "On the generality of refined algebraic quantization," 
Class. Quant. Grav. 16, 2479 (1999) |arXiv:gr-qc/9812024] . 

[59] D. Giulini and D. Marolf, "A Uniqueness theorem for constraint quantization," 
Class. Quant. Grav. 16, 2489 (1999) ,arXiv:gr-qc/9902045,. 

[60] D. Giulini, "Group averaging and refined algebraic quantization," Nucl. Phys. Proc. 
Suppl. 88, 385 (2000) |arXiv:gr-qc/0003040| . 

[61] T. Thiemann, "Quantum spin dynamics. VIII. The Master constraint," Class. 
Quant. Grav. 23, 2249 (2006) |arXiv:gr-qc/0510011 . 

[62] M. Han and T. Thiemann, "On the Relation between Rigging Inner Product and 
Master Constraint Direct Integral Decomposition," J. Math. Phys. 51, 092501 
(2010) ■;arXiv:0911.343r [gr-qc]]. 

[63] A. Laddha and M. Varadarajan, "The DifFeomorphism Constraint Operator in Loop 
Quantum Gravity," Class. Quant. Grav. 28, 195010 (2011) |arXiy:1105.0636 [gr-qc]]. 

[64] D. Marolf, J. Mourao and T. Thiemann, "The Status of diffeomorphism su- 
perselection in Euclidean (2-1-1) gravity," J. Math. Phys. 38, 4730 (1997) 
^arXiv:gr-qc/9701068 . 

[65] H. Sahlmann, "Exploring the diffeomorphism invariant Hilbert space of a scalar 
field," Class. Quant. Grav. 24, 4601 (2007) [arXiv:gr-qc/0609032| . 

[66] W. Fairbairn and C. Rovelli, "Separable Hilbert space in loop quantum gravity," J. 
Math. Phys. 45, 2802 (2004) |arXiv:gr-qc/0403047 . 

[67] C. Fleischhack, "Representations of the Weyl algebra in quantum geometry," Com- 
mun. Math. Phys. 285, 67 (2009) |, arXiv:math-ph/0407006] . 

[68] M. Varadarajan, "Towards new background independent representations for loop 
quantum gravity," Class. Quant. Grav. 25, 105011 (2008) [arXiv:0709.1680 [gr-qc]]. 

[69] M. Dziendzikowski and A. Okolow, "New diffeomorphism invariant states on a 
holonomy-fiux algebra," Class. Quant. Grav. 27, 225005 (2010) |arXiv:0912.T278l 

[math-ph]]. 

[70] T. A. Koslowski, "Dynamical Quantum Geometry (DQG Programme)," 
LarXiy;0709.3465 [gr-qc]. 

[71] H. Sahlmann, "On loop quantum gravity kinematics with non-degenerate spatial 
Class. Quant. Grav. 27, 225007 (2010) |arXiv:10 06.0388 [gr-qc]]. 

[72] T. Koslowski and H. Sahlmann, "Loop quantum gravity vacuum with nondegenerate 
geometry," arXiv:11 09.4688 [gr-qc]. 

[73] T. Thiemann, "Quan tum spin dynamics (QSD)," Class. Quant. Grav. 15, 839 (1998) 
[arXiv:gr-qc/ 9606089 . 

[74] T. Thiemann, "Quantum spin dynamics (qsd). 2," Class. Quant. Grav. 15, 875 
(1998) ,arXiv:gr-qc/9606 090 . 

[75] T. Thiemann, "QSD 3: Quantum constraint algebra and physical scalar prod- 
uct in quantum general relativity," Class. Quant. Grav. 15, 1207 (1998) 
[arXiv:gr-qc/ 97050l7 . 

[76] M. P. Reisenberger and C. Rovelli, "'Sum over surfaces' form of loop quantum 
gravity," Phys. Rev. D 56, 3490 (1997) [arXiv:gr-qc /9612035. 

[77] C. Rovelli, "Quantum gravity as a 'sum over surfaces'," Nucl. Phys. Proc. Suppl. 
57, 28 (1997). 

54 



[78] J. Lewandowski and D. Marolf, "Loop constraints: A Habitat and their algebra," 
Int. J. Mod. Phys. D 7, 299 (1998) |arXiv:gr-qc/9710016| . 

[79] A. Laddha and M. Varadarajan, "Polymer quantization of the free scalar field and 
its classical limit," Class. Quant. Grav. 27, 175010 (2010) arXiv:1001.3505 [gr-qc]]. 

[80] A. Laddha and M. Varadarajan, "The Hamiltonian constraint in Polymer 
Parametrized Field Theory," Phys. Rev. D 83, 025019 (2011) [arXiv: 1011 .24631 [gr- 

qc]]. 

[81] A. Perez, "On the regularization ambiguities in loop quantum gravity," Phys. Rev. 
D 73, 044007 (2006) | arXiv:gr-qc/050911^ |. 

[82] T. Thiemann, "Reduced phase space quantization and Dirac observables," Class. 
Quant. Grav. 23, 1163 (2006) [arXiv:gr-qc/0411031 . 

[83] A. S. Vytheeswaran, "Gauge unfixing in second class constrained systems," Annals 
Phys. 236, 297 (1994). 

[84] K. Giesel and T. Thiemann, "Algebraic Quantum Gravity (AQG). II. Semiclassical 
Analysis," Class. Quant. Grav. 24 (2007) 2499 gr-qc/0607100|. 

[85] M. Domagala, K. Giesel, W. Kaminski and J. Lewandowski, "Gravity quantized: 
Loop Quantum Gravity with a Scalar Field," Phys. Rev. D 82 (2010) 104038 
[arXiv: 1009.2445 [gr-qc]]. 

[86] J. Brunnemann and T. A. Koslowski, "Symmetry Reduction of Loop Quantum 
Gravity," Class. Quant. Grav. 28, 245014 (2011) |arXiv:1012.005^ l [gr-qc]]. 

[87] C. Fleischhack, "Loop Quantization and Symmetry: Configuration Spaces," 
[irXiviipiO.0449 [math-ph]. 

[88] T. Thiemann, "Gauge field theory coherent states (GCS): 1. General properties," 
Class. Quant. Grav. 18, 2025 (2001) [arXiv:hep-th/0005233| . 

[89] T. Thiemann, "Complexifier coherent states for quantum general relativity," Class. 
Quant. Grav. 23, 2063 (2006) [arXiv:gr-qc/020603r. 

[90] B. Bahr and T. Thiemann, "Gauge-invariant coherent states for loop quantum 
gravity. II. Non-Abehan gauge groups," Class. Quant. Grav. 26, 045012 (2009) 
[arXiv:0709.4636l [gr-qc]]. 

[91] C. Flori and T. Thiemann, "Semiclassical analysis of the Loop Quantum Gravity 
volume operator. I. Flux Coherent States," arXiv:0812.1537 [gr-qc]. 

[92] C. Flori, "Semiclassical analysis of the Loop Quantum Gravity volume operator: 
Area Coherent States," |arXii:0904. 1303 [gr-qc]. 

[93] L. Freidel and E. R. Livine, "U(N) Coherent States for Loop Quantum Gravity," J. 
Math. Phys. 52, 052502 (2011) ar Xiv:1005 .2090 [gr-qc]]. 

[94] F. Girelli and E. R. Livine, "Reconstructing quantum geometry from quantum in- 
formation: Spin networks as harmonic oscillators," Class. Quant. Grav. 22, 3295 
(2005) arXiv:gr-qc/0501075 . 

[95] L. Freidel and E. R. Livine, "The Fine Structure of SU(2) Intertwiners from U(N) 
Representations," J. Math. Phys. 51, 082502 (2010) [arXiv:091 1.3553 [gr-qc]]. 

[96] L. Freidel and S. Speziale, "Twisted geometries: A geometric parametrisation of 
SU(2) phase space," Phys. Rev. D 82, 084040 (2010) [arXiSOOl.2748 [gr-qc]]. 

[97] C. Rovelli and S. Speziale, "On the geometry of loop quantum gravity on a graph," 
Phys. Rev. D 82, 044018 (2010) |arXiv:1005.2927l [gr-qc]]. 

[98] E. Bianchi, P. Dona and S. Speziale, "Polyhedra in loop quantum gravity," Phys. 
Rev. D 83, 044035 (2011) [arXiv: 1009.3402 [gr-qc]]. 

[99] E. Bianchi and H. M. Haggard, "Discreteness of the volume of space from Bohr- 
Sommerfeld quantization," Phys. Rev. Lett. 107, 011301 (2011) [arXiv: 1102. 54391 

[gr-qc]]. 

55 



[100] E. F. Borja, L. Freidel, I. Garay and E. R. Livine, "U(N) tools for Loop Quan- 
tum Gravity: The Return of the Spinor," Class. Quant. Grav. 28, 055005 (2011) 
arXiv:1010.5451 [gr-qc]]. 

[101] E. R. Livine and J. Tanibornino, "Spinor Representation for Loop Quantum Grav- 
ity," J. Math. Phys. 53, 012503 (2012) |arXiv:1105.3385 [gr-qc]]. 

[102] E. R. Livine, S. Speziale and J. Tambornino, "Twistor Networks and Covariant 
Twisted Geometries," Phys. Rev. D 85, 064002 (2012) [arXi v:1108.0369i [gr-qc]]. 

[103] E. R. Livine and J. Tambornino, "Loop gravity in terms of spinors," larXiv:1109.3572l 

[gr-qc]. 



56 



